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ABSTRACT 


rAnrnn  iit  nr 


a  ruKiKAiN  iv  program  has  been  written  which  uses  an  iterative  tech- 
nique  to  model  ion  transport  and  chemistry  in  a  drift  tube.  The  program, 
named  MULTIGATE,  uses  finite  difference  calculus  to  generate  time  arrival 

etAri  Ae  of  c  1  r\ni  nf  c  i  n  f  Via  /It1!  *Pt  cnaro  for  tlio  PQ  CP  n  tllTPP  !  OT! 

uxg  VV4i  XV  J  u.  W  V  V  V  X  U>  X  VX  II  XII  V1IV  V4X  XX  W  V  xwx  V**  V  v  w*  w 

species  undergoing  populating  and  depopulating  reactions  and  at  the  same 
time  drifting  with  the  electric  field  and  diffusing  in  the  radial  and 
dirertinrn;.  Tn  facilitate  comnarison  with  drift  tube  data,  the 


computer  output  mciuaes  plots  or  rne  ionic  currents  Versus  time. 
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GLOSSARY 


Throughout  the  glossary  the  symbol  A  is  used  to  refer  to  ion  species 
A.  The  symbols  B  or  C  may  be  substituted  to  form  the  FORTRAN  variable 


n allies  applying  to  ion  species  B  gt  C . 


Til/  aiji  c  a 


♦  etnmKrtl  TT1/H  mov  V\  c 

uuw  rv  i -/v 


replaced  by  the  symbol  "1,"  "2,"  "3,"  "4,"  or  a  blank  to  refer  to  the 
four  intermediate  and  the  final  gates  in  the  list  of  FORTRAN  variable 


nomac 


A 

AP 

n  a 

rn 

APP 

AGAMk 

DA 

n&D 

um\ 

DAZ 

DNA 

DR 

rw 

UU 

DT 

E 

EON 

IAk 

INGAk 

Kk 

TV 


MAAAK 


calculated  value  of  the  density  at  time  T  of  ion  species  A 
time  derivative  of  the  density  at  time  T  of  ion  species  A 

A  J  a  ^  1  ••  A  A  IT  A  J  AM  A  4  A  ^  4"  4  tn  +S.  Tx  IYT  A  4*  1  P  1A  A/*  i  AP  A 

piCUlL  L  cu  value  U1  LHC  UCUDll/  ac  UJLllIt;  1TU1  Ui  lull  rx 

time  derivative  of  the  predicted  density  at  time  T+DT  of  ion 

species  A 

array  of  length  MAXZ  in  which  the  currents  due  to  ion  species  A 
at  gate  k  are  stored 

total  distance  that  the  cell  space  for  species  A  has  drifted 

ol  /Ti  ffnei  rm  rnoffi  ri  out  fnr  inn  cnnninc  A 

J.  UVIXU1  VIXi.  X  XVII  WVX  X  X  VX  W  *  wx  x  Vi*  wx  w 

axial  diffusion  coefficient  for  ion  species  A 

initial  density  of  species  A  measured  at  the  peak  of  the 
distribution 

grid  size  in  the  radial  direction 

/w»  -i  a  i  -r  a  -i  a  f  Vi  a  o  vi  o  1  ■?  ta  An 
^11U  J1  111  LUU  OAiai  U  J.X  X  VSI1 

time  step 

electric  field  value  (volts/cm) 

2 

ratio  of  electric  field  to  gas  number  density  (volts  cm  ) 

counter  for  storing  the  currents  at  gate  k  for  ion  species  A 

switch  which  indicates  when  distribution  A  has  passed  gate  k 

represents  the  distance  from  leftmost  cell  of  the  cell  space  A  to 
gate  k  (cell  space  assumed  drifting  from  left  to  right) 

represents  the  distance  from  the  leftmost  cell  of  the  cell  space 
C  to  gate  k 

nAn-nne  nn  +  e  4*  1a  a  A  -i  etonon  ft*nm  tV>o  1  a-pfmnct*  col  1  nf  tKo  AO  1  1  cno  AO 

X  upl  UOUU  1-0  L11U  UXJtUUWV  l  X  will  Vliv  xvx  WUWJ  v  WAX  wx  waxvr  w  x  x 

B  to  gate  k 

internal  switch  which  indicates  when  the  cell  space  for  species 
A  is  traversing  gate  k 


Q 


GLOSSARY  (Contd) 


MAXZ  number  of  cells  in  the  axial  direction 

PI  pressure  of  the  reactant  gas  through  which  the  ions  are  drifting 

(Torr) 

P 2  pressure  of  the  buffer  gas  through  which  the  ions  are  drifting 

(Torr) 

RAB  frequency  of  the  reaction  producing  ion  species  A  from  ion 
species  B 

RO  Gaussian  halfwidth  in  radial  direction  of  initial  ion  species 

distributions 

SING  indicates  whether  rescaling  is  necessary  for  species  A  in  the 
radial  or  axial  direction 

T  total  drift  time 

TAk  array  of  length  MAXZ  containing  corresponding  times  for  AGAMk 

TEMP  temperature  of  the  system 

TING  indicates  whether  rescaling  is  necessary  for  species  B  in  the 
radial  or  axial  direction 

VA  drift  velocity  of  ion  species  A 

XKA  mobility  of  ion  species  A 

XKOA  reduced  mobility  of  ion  species  A 

XLk  position  of  gate  k 

XN1  density  of  gas  1 

XN2  density  of  gas  2 

YKAB  rate  constant  for  ion  species  A  converting  to  B 

ZING  indicates  whether  rescaling  is  necessary  for  ion  species  C  in  the 
radial  or  axial  direction 

ZO  full  Gaussian  width  in  the  axial  direction  of  the  inital  ion 

species  distribution 
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I.  INTRODUCTION 


The  performance  of  communications  systems,  radars,  and  other  defense 
systems  are  often  predicted  using  phenomenology  codes,  some  of  which  employ 
simplified  chemistry.  Extensive  chemical  codes  of  the  lower  ionosphere 
are  being  used  to  check  the  reliability  of  the  simplified  chemistry.  The 
development  of  these  large  chemical  codes  requires  both  in  situ  measure¬ 
ments  of  ionic-  and  neutral-constituent  number  densities  and  laboratory 

m  non  tVAm  a  rt'T  +-U  /\  1  /%irnr>4-  via  n  a  vi  a  +  a  r>  Hn  a  /~\  £  +  U  a  1 

iii&aoui  o  \j  x.  luc  i  cicvaiiL  icaLLiuu  i  auc  Luiioiamo  .  uuc  ui  luc  xclukjxcl  — 

tory  tools  that  is  particularly  well  suited  for  studying  ion-molecule 
reactions  is  the  drift  tube.  One  of  these  is  now  in  operation  at  the 

Rfl  1  1  i  r  R  a  t  r  h  I.ahnratnri  pq  fRRT.'l 


A  typical  drift  tube  configuration  is  shown  in  Figure  1.  Ions  are 
produced  in  the  ion  source.  An  axial  electric  field  causes  them  to  drift 
down  the  tube,  and  most  of  them  impact  the  end  wail.  Some  of  the  ions 
are  swept  through  the  exit  aperture  into  the  mass  spectrometer.  Typically 
the  source  runs  continuously,  with  one  of  the  pairs  of  electrical  grids, 
called  "shutters,"  serving  as  the  point  of  origin  for  a  swarm  of  ions. 

That  shutter,  the  "starting  shutter,"  is  electrically  biased  to  stop  the 
drift  of  ions;  then  it  is  opened  by  applying  a  pulse  of  dc  voltage  for  a 

■Pout  mi  r'rncnpATi/lc  t  a  oil  Ai.r  o  c«o  +  i  o  1  1  x;  email  nmnn  P  i  nr>c  t  a  ctort  A  Aum 
i.  vn  v^v/wiiuj  lu  ax  xv/n  a  jpu  c.  a.  c*  x.  x  j  umux  x  giva^  w  x  xaiu  cw  o  xui  c  avnit 

the  tube.  As  the  cloud  of  ions  drifts,  it  diffuses  both  axially  and 
radially,  and  the  ions  react  with  the  gas  molecules  in  the  tube  to  form 
other  ions  which  have  their  own  drift  velocities  and  diffusion  coeffi¬ 
cients.  The  mass  spectrometer  is  used  to  select  ions  of  one  mass  number 

from  the  ions  that  emerge  from  the  exit  aperture,  and  the  signals  due  to 

the  arrival  of  ions  of  that  mass  are  sorted  and  stored  in  a  multichannel 

nnn  1  \/7or  o  c  a  Pi  mof  i  a«  aP  a  ti  mo  ci  nno  t  V\o  c  +  ort  i  nrr  c  Vuif  +  or>  uo  c  tm  i  1  co^ 

ana  x  j  x  a  j  a  x  cxi  l  v_  cxvn  v/x  cue  cxiue  oxuvv  cue  j  cux  c  x  ug,  uuuccex  nuo  pux  jeu 

on.  Typically,  fewer  than  one  ion  per  ion  group  negotiates  the  drift 
tube  and  mass  analysis  chamber  and  is  counted,  so  that  the  starting 
shutter  must  be  pulsed  many  times  before  statistically  significant  arrival 
time  profiles  for  the  several  ions  are  accumulated.  A  complete  set  of 
these  profiles  for  two  gas  pressures  at  a  convenient  value  of  E/N,  where 
E  is  the  drift  field  strength  and  N  is  the  gas  number  density,  contains 
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DRIFT  TUBE 


Figure  1.  Typical  Drift  Tube  Configuration 


the  information  necessary  to  deduce  the  drift  velocities  (and  thus 


mr»V\i  litiocl  *yr\A  1-  Vi  o  roontinn  t*q+,o  rrmctantc 

1IIW  UXAXWXV«J  ^  WIIU  V  1  1  V  i.  VU  V  W  O.V  II  A  W  Wl&«^  V  VklA  V  ( 


llnf  nrtimQf  a  1  vr 

WAiA.  V/  W  V»1IV»  V  W 


constants  and  some  of  the  drift  velocities  cannot  be  easily  deduced  from 
the  arrival  time  curves.  One  approach  is  to  compare  measured  arrival 
time  curves  with  curves  calculated  for  a  variety  of  Darameters  to  insure 

___  _  _  __  __  - -  ___  __  /  i 


that  the  uCSIi  Cd  physical  Quantities  ax  e  UiilQuc  1/  determined .  In  gen- 

erating  the  calculated  arrival  time  curves  one  must  take  into  account  all 
the  known  physical  phenomena  which  will  affect  the  shapes  of  the  curves, 

i  r\nr  oil  norci  kl  a  v  ao  i  at>c  oil  ava  i  nnr  i  r\r\  c  A  onmni  i  +  av*  nrnrrr  om 

i.Jl^iUUXUg  ux  x  i/x  v  x  uuw  ux  vx.  ux  X  VlilVX  ^xu^  XVU«7  •  n  Wlll^/U  WVX  UJIA  y 

MULTIGATE,  has  been  developed  at  the  BRL  which  permits  numerical  modeling 
of  the  drift,  diffusion,  and  reactions  of  several  ions  in  a  drift  tube. 


MULTIGATE  is  a  FORTRAN  IV  program,  originally  conceived  of  by  Kregel, 
that  uses  an  iterative  numerical  technique  for  drift  tube  modeling. 
MULTIGATE  uses  finite  difference  calculus  to  generate  time  arrival  his¬ 


tories  at  several  points,  or  •gates,"  m  tne  amt  space  tor  several  ionic 
species.  The  modified  Euler  method  is  employed  for  performing  the  actual 
numerical  integration  phase  of  the  program.  At  present,  the  program  is 


U/m  i  /»  U  1-i  m  •!  4> 

UlUU^ll  UUL  llllliigu  LU 


hree  ion  spec 1 0 s  and  four  inter- 


mediate  gates,  and  one  end  gate,  hence  the  name  MULTIGATE. 


In  the  program  finite  cylindrical  geometry  (r,z)  is  assumed  with 
the  drift  space  divided  into  four  concentric  shells  of  equal  thichness  DR 
which  are  each  in  turn  subdivided  along  Z  into  sections  of  length  DZ. 

Ion  densities  are  defined  at  the  center  of  each  of  the  resulting  Mringsff 
and  are  assumed  to  vary  linearly  between  the  centers.  Each  ring  defines 
a  cell  for  computational  purposes.  It  has  been  found  that  63  divisions 
along  the  z  axis  (making  a  total  of  252  cells  in  all  for  each  species')  are 


adequate  to  effect  a  reasonable  compromise  between  truncation  error  and 
computational  time.  The  number  of  cells  in  the  cell  space  for  each  ion 
is  arbitrary  provided  that  it  is  sufficiently  large.  One  unique  feature 
of  the  program  which  allows  a  significant  reduction  in  the  number  of  cells 
required  is  to  have  the  ceil  space  for  each  ion  "drift"  with  the  drift 
velocity  of  the  corresponding  ion  species.  This  permits  the  peak  of  each 
ion  cloud  to  remain  more  or  less  centered  in  its  corresponding  cell  space. 
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MAIN  PROGRAM 


II . 


The  primary  purpose  of  the  main  program  is  to  control  the  flow  of 

rnrrmi 1 1 a t  i  rm  anH  tn  nnntrnl  all  niitnnt  .  Tt  initial  lv  the  limitina 

- L - -  —  —  — * - r - -  —  - j  — -  “  -  - o 

values  of  the  variables,  determines  the  size  of  each  cell,  sets  the  gate 
distances,  and  generates  an  initial  distribution  for  each  ion  species  in 
the  aDDroDriate  cell  soaces.  At  the  beginning  of  each  run.  the  following 

XXX  X  ^  ^  ^ 

parameters  are  initial! zed ; 


Symbol  Used 
in  the  Program 


Symbol  Definitions 


XKOA,  XKOB,  XKOC 


The  reduced  mobilities  of  ion  species  A,  B, 
and  C,  respectively  (cm^V^sec^)  . 


hma  nMR  r.wr 


Initial  densities  of  species  A ,  B,  and  C 
measured  at  the  neak  of  the  distribution. 


ZO,  RO,  XL 


Full  Gaussian  width  of  the  initial  ion 
distributions  in  the  axial  direction,  the 
Gaussian  half  width  of  the  initial  ion 
distributions  in  the  radial  direction,  and 
total  length  of  drift  space  (initially,  the 
ion  distributions  are  centered  at  Z=0  and 
later  drift  to  Z=XL) . 


PI,  P2 ,  TEMP,  EON 


Pressure  of  the  reactant  gas  through  which 
the  ions  are  drifting  (Torr) ,  pressure  of 


the  buffer  gas ,  temperature  of  the  gases 


rm 


( o  c  oniio  1  ^  o-nrl  rot  i  n  r\4- 

V\^UUX  J  ,  UUU  1.  J.V/  Vi- 


field  to  total  pas  number  densitv  fVcnrl . 


> 1 ortri r 


YKAB 


Rate  constant  for  reaction  which  forms  B 
from  A,  or  "A  B. " 


\nrr\  a  \m^  a  /n  \M/n  a 

I NDA ,  INAt,  INCA, 

YKBC ,  YKCB 


Rate  constants  fur  reactions  B  A,  A  ^  C , 
C  A,  B  C,  and  C  -t  B. 


The  main  program  also  sets  certain  internal  switches  and  counters 

lir  a  J  i  4-  Vi  l  f  n  -v»  -1  /Mir  riiV\Y»/Mi  +  i  n  A  r*  4-  rt  /I  of  n  vm  1  n  Q  4  K  n  C"  4  T  f  Q  nf  4  V>  O  r*  riTTI  Y~\  1  l4  Q  4  1  r\n 

UJCU  JL11  Lilt;  V  CXL  1UUO  OUU1  UULXIIUU  IV  UV  LV^imi.HU  UllV  J  LU-UV  \j  x.  CilW  WWIH^U  v.  c*  v-  x.  ^41 

(defined  below  for  gate  1) : 


MA  Y  A  1  y  A  V  D  1  \,1  A  V  P  1  T  nf/ivno  1  ct.Ti  4  c>Vi  oc  ■Pnr'  i  n/1  i  ^  o  4  t  r\  rr  wrVi  cm  4  Vl  Ci 

lTLrv/vrvx  ,  ,  I'IAAVj  X  xnu/Uiax  J  i  l  x.  \j  x  x.  nu  x.  v-a  u  j-  ug  nuvu  tuv 


cell  snace  for  snecies  A.  B.  and  C.  resnec- 

- —I - - -X -  -  *  -  ’  -  '  1 

tively,  is  traversing  gate  1. 


IAA1,  IBB1 ,  ICC1  Counters  used  to  indicate  where  the  values 

for  the  currents  and  corresponding  times 
for  ion  species  A,  B,  and  C,  respectively, 
are  stored. 


INGA1,  INGB1,  INGC1  Internal  switches  for  indicating  when  the  cell 

spaces  for  species  A,  B,  and  C,  respectively, 
have  traversed  gate  1. 

Kl,  Ml,  LI  Counters  used  to  indicate  the  position  (in 

units  of  cell  thickness)  of  gate  1  with 
respect  to  the  cell  spaces  for  species  A, 

B,  and  C,  respectively. 

The  program  contains  similar  internal  switches  and  counters  for  the 
last  four  gates. 


III.  SELECTION  OF  DT 


Throughout  the  computation  the  time  step,  DT,  used  in  the  integration 
of  the  ion  densities  is  maintained  at  each  time  step  at  a  value  that  will 
assure  stability  and  minimize  run  time. 

The  maximum  time  step  permitted  by  diffusion  for  each  of  the  three  ion 
species  is  computed  by  the  following  method.  The  particle  current  density 
J  due  to  diffusion  is  defined  vectorially  as 

J  =  -DVn  (1) 

where  D  is  the  diffusion  coefficient  and  n  is  the  particle  number  density. 
Since  cylindrical  symmetry  is  assumed  in  our  case,  J  may  be  approximated 
using  finite  difference  calculus  as 


-r  „  An  —  an  — 

J  as  -D  —  1  -  D  - —  1 

z  Az  z  r  Ar  r 


«  J  i  +  J  i 

z  z  r  r 


(2) 

(3) 


in  which  i^  and  i^  are  unit  vectors  in  the  z  and  r  direction,  respectively, 

and  An,  Ar,  and  Az  are  finite  differences  denoting  differences  in  density, 

in  distance  in  the  i  direction  and  in  distance  in  the  i  direction, 

r  z  ’ 

respectively. 

In  the  following  discussion  we  shall  consider  a  cell  of  thickness 
Ar  and  length  Az  located  at  rQ  and  zq.  From  the  continuity  equation 
(also  cast  into  a  finite  difference  equation)  we  know  that  the  particle 
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current  density  evaluated  at  a  particular  surface  of  a  cell  times  the 
surface  area  summed  over  all  cell  surfaces  when  divided  by  the  cell 

irrtliimrt  i.t  A  1  1  <vrtr\-»*AV  i  mo  +  a  +U  a  nArrn  +  -lira  rvf  +■  l  m  /-v  vo  +  a  r\  4-  onrra  /~\  -P  f  V>  a 

v  uiuiuu  m  x  i  i  app  l  UAxma  Xiiv  iibgotivt  ux  xuv^  txmt  x  a  tx  w  x_  xiiaugx  wx.  xnx 

density,  n,  in  that  cell.  Current  is  defined  to  be  positive  if  it  is 

directed  in  the  direction  of  the  appropriate  unit  vector.  Thus  the 

instantaneous  time  dpri vat i vp  of  n  rnav  hp  writtpn  a<; 


-A(zJ  Jjz J  +  A(z„  +  Az)  JJz„  +  Az)  -  A(rJ  (rj  +  A(r^  +  Ar)  J^(r^+Ar) 


'dt 


2Trr  ArAz 
o 


(4) 


where  A(z  )  and  J  (z  )  denote  the  cell  area  defined  by  the  surface  for 
o  z  o  1 

which  z=z,r  <  r  <  r  +Ar  and  the  current  through  that  surface, 

0  0  0 

respectively,  and  where  the  other  areas  and  currents  are  defined  in  a 
similar  manner.  All  areas  are  assumed  to  be  positive  in  extent. 

We  shall  let  An  denote  the  amount  by  which  the  ion  number  density 
in  the  ceil  differs  from  the  "average”  density  in  each  of  the  surrounding 
cells,  and  shall  assume  that 


fin 


J„(zJ  «  -JJz„  +  Az)  =  -D„  i^- 


and 


J  (r  )  -  -J  (r  +  Ar)  h  -D 


Substitution  yields 


an 

J  X 
U  L 


Afz  1  +  Afz  +  Azl 

D 

Afr  1  +  Afr  +  Arl 

'  O'  '  0  'J 

a  z  r 

'  0"  "0  J 

(5) 


2Tir  ArAz 
o 


If 


A(zq)  «  A(zq  +  Az)  *=  2Trr  Ar 


then 


j  D  An  D  An 

an  ^  z  ^  r 


i6j 


c - _ r  £  \ _ i_ _ _  j: 


l.  i _ r* _ j 


LjiL^Lid nun  ^oj  empiiasxzes  Lne  iaci  tnat  ir  an  is  positive,  xne  xime  rare  ox 
change  of  the  ion  number  density  is  negative;  that  is,  there  is  a  net  loss 
of  ions  from  the  cell.  If  An  is  assumed  to  vary  linearly  in  time  in  the 

1  n  4-  AW  r  n  1  A  4-  cm/1  -l  ^  /»/M  m  1  -5  «/•  /-N  4»  a  n  J  «■»  r.  4  ^  4  /N  n  M  ^  L  ~  ^1.  ~  —  ~  ~  1  1  —  1  _ 

-luc^a  v  a.  J.  <->  L.  auu  ii  UCtWCCIl  L11C  UCIWlllCi  XII  LI1C  ULI1CX  tCiib  ib 

neglected,  the  time  increment  associated  with  a  change  of  An  in  the  den¬ 
sity  in  a  cell  can  be  expressed,  in  terms  of  finite  differences,  as 


At  « 


.  /D 


D  \ 

..  .  A-) 

^\Ar"  Az  7 


-1 


(7) 


Limiting  At  then  limits  the  possible  An  in  anv  one  time  step.  To  insure 

x  #  j. 

C  t  O  V\  -ill  /~\  ■C  V»  /\  r>  rt  1  i  i  ♦  -I  A  n  A  4-  •!  <•  n  n4-  4>  A  A  7C  ^  n  4-  L  1  n  ^  1  .  .  A  1  _  J  »  ^  4. 

J  I~c xu/xxxu/  XllC  DU1UU.U11,  LiL  ib  bCL  IU  U.  /J  XXHICS  LIIX25  VaXUC  XJ1  UiUCI  UU 

damp  out  any  Mnoise"  on  n. 

AO  1  A1.1  n4>4  4-  « 1,  ~  A.  r»  ^  «  J  -!  „  4 - V  T  ^  ^ 

\juiigiit  uaiLUiaiiuiib  lciin.c  piatc  clz>  ^ciia  uun  a  uibiautc  Aa  xn  uic 

x-direction,  that  is,  as  they  pass  the  imaginary  end  plate.  Consideration 
must  be  made  of  the  effects  of  the  distance  that  each  cell  space  drifts 


during  a  single  time  step.  If  the  time  step  were  so  large  th  — t*  -  p«>i  i 


w-  V  V  1.  X 


i  •  r*  . 


space  onrteo  a  distance  equal  to  several  ul ,  mere  would  oe  tewer 
current  calculations  than  there  were  axial  cells  for  that  ion,  thus 
diminishing  the  desired  output  of  the  program.  Therefore  we  calculate  the 

m  n  v  i  mi  iwi  r*  i  rt  ^  ^  ^  1%  ^  ^  4  ^  4-  «%  ^  ^  J  ^  J  i? x.  %  .  A  1  i  x.  _  _  ^  r  ^  1  _  _  • 

jixcxw  xiiiciiu  wx  cue;  Lxiuc  otep  penult  ucu  u/  cue;  unit  vci  uci  ties  or  me  ions 

by  the  formula 


rr  V  n7 
t  *■'  L/trf 

At  =  - -  , 

v 

max 


(8) 


where  v  is  the  drift  velocity  of  the  fastest  ±un,  g  is  a  control  factor 
max 

equal  to  the  number  of  cells  which  may  drift  past  the  end  plate  in  one  time 
step,  and  DZ  is  the  axial  grid  size.  This  procedure  insures  that  there  will 


Ka  O  Kai  1  f  MA  Y7  /  rr  All  w  rtr»+-  ao  1  / 
u'*  UUUUL  riry.j\i*/  g  cunciii  can, 


.  . . - „  1  1  .  - 4. 

nc  UbUdiiy  bCL 


g  =  1.5. 
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Associated  with  each  reaction  is  a  rate  constant,  e.g.,  YKAB,  and  a 
frequency  of  reaction,  e.g.,  RAB.  The  quantity  1.0/RAB  is  the  time 
constant  for  the  reaction  whose  rate  constant  is  YKAB.  The  At  for  this 
reaction  is  set  at  0.8  x  1.0/RAB,  again  to  insure  stability  of  the 
solution. 

The  program  sets  DT  equal  to  the  minimum  time  step  calculated  for 
diffusion,  drift,  or  any  of  the  reactions. 

IV.  INITIAL  DISTRIBUTION 

Initially  the  program  assumes  Gaussian  distributions  m  both  the 
radial  and  axial  directions  for  the  densities  of  all  of  the  ion  species. 
The  distribution  of  ion  species  A  is  calculated  from  the  relationship: 

A(r ,  z)  =  DNA  x  EXp|-ar2  -  S(z  -  zq)  2J  ,  (9) 

where  DNA  is  the  initial  value  of  A(r,z)  at  the  peak  of  the  distribution, 

a  is  calculated  by  the  program  from  RO,  the  Gaussian  halfwidth  in  the 
radial  direction,  to  make  the  density  in  the  fourth  radial  zone  approxi¬ 
mately  1/e  its  value  at  the  center,  and  B  is  calculated  from  ZO,  the  full 

Gaussian  width  in  the  axial  direction,  to  make  the  full  width  of  the 

density  profile  in  the  longitudinal  direction  approximately  ten  cells. 

The  parameter  ZO  insures  that  the  initial  distribution  is  centered  in  the 
cell  space.  The  program  then  initializes  the  distributions  for  ion 
species  B  and  C  according  to 


B(r, z) 

=  A(r , z)  x  DNB/DNA 

(10) 

C  (r ,  z) 

=  A(r, z)  x  DNC/DNA  , 

(11) 

so  that  all  species  have  initially  the  same  distribution  shapes. 

V.  DERIVATIVE  SUBROUTINE  -  S 

Subroutine  S  is  used  to  calculate  time  derivatives  of  the  ion  number 
densities  and  of  predicted  ion  number  densities  in  the  cells  in  the  three 
inner  radial  shells.  Effects  of  drift,  diffusion,  and  chemistry  are 
handled  in  S. 
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The  calculation  begins  by  approximating  the  current  density  at  the 
four  boundaries  of  each  cell  using  second  order  difference  calculus. 

These  four  current  densities  are  then  multiplied  by  the  corresponding 
incremental  areas  and  summed.  This  resulting  sum  is  then  divided  by  the 
incremental  volume  of  the  corresponding  cell  to  yield  the  time  rate  of  change 
of  density  due  to  diffusion.  Effects  of  mobility  are  taken  into  account 
automatically  since  the  cells  are  allowed  to  "drift"  with  the  drift  velocity 
of  the  corresponding  ionic  species. 

The  time  rate  of  change  of  the  density  due  to  chemistry  for  each  ionic 
species  at  each  cell  is  easily  calculated  by  multiplying  the  rate  constants 
by  the  appropriate  densities  of  the  interacting  ionic  species  evaluated  at 
the  center  of  that  cell.  First  order  extrapolation  techniques  are  used 
for  this. 


Thus,  in  generating  the  derivatives  due  to  diffusion  the  subroutine 
calculates  the  current  densities  times  the  incremental  area  in  consecutive 
cells  in  the  radial  and  axial  directions.  The  following  names  are  used 
in  subroutine  S: 


HA's  Represent  densities  of  A. 

GA's  Denote  product  of  current  density  and  area  or  the 

number  of  ions/sec  passing  through  an  elementary  area. 

HAP  Density  of  A  in  cell  (i=l,2.3.  or  4:  k=1.2.3 .  or 

MAXZ) . " 

HA  Density  of  A  in  cell  (j,k+l),  in  increasing  axial 

direction. 

HAR  Density  of  A  in  cell  (j+l,k),  in  increasing  radial 

direction. 

GAZP  Current  flowing  into  cell  (j  ,k)  in  axial  direction. 

GAZ  Current  flowing  out  of  cell  (j  ,k)  in  axial  direction. 

GARP  Current  flowing  into  cell  (j ,k)  in  radial  direction. 

GAR  Current  flowing  out  of  cell  (j ,k)  in  radial  direction. 


In  addition  to  the  ion  diffusion,  the  program  must  also  calculate 
the  effects  of  chemical  reactions  upon  the  three  ions.  Since  the  cell 
spaces  of  the  three  ions,  in  general,  do  not  drift  at  the  same  speed,  the 
program  must  keep  track  of  the  relative  shift  of  the  cell  spaces  in  order 
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to  "bring  together"  the  points  of  equal  drift  distance  of  the  three  ions. 
To  accomplish  this,  the  main  program  records  the  drift  distances  of  the 
three  cell  spaces  in  DA,  DB,  and  DC.  Each  drift  distance  is,  of  course, 
the  sum  over  time  of  the  corresponding  drift  velocity  times  the  appro¬ 
priate  time  step. 

VI.  SCALING  -  BANGR,  RAP,  BANGZ,  AND  ZAP 

During  the  calculation,  tests  are  made  each  time  step  to  determine 
whether  rescaling  is  necessary  in  either  the  radial  or  axial  direction. 
The  subroutine  BANGR  is  used  to  determine  whether  rescaling  is  necessary 
in  the  radial  direction.  Since  the  ion  distributions  tend  to  remain 
centered  in  their  respective  cell  spaces,  it  is  only  necessary  to  test 
in  a  plane  perpendicular  to  the  axis  halfway  down  the  cell  space  of  each 
ion.  Rescaling  in  the  radial  direction  occurs  only  when,  for  any  ion 
species,  the  ratio  of  the  density  in  the  cell  in  the  outer  radial  shell 
to  the  cell  in  the  innermost  shell  is  greater  than  0.9,  a  number  chosen 
by  trial  and  error.  Number  densities  corresponding  to  an  expanded  radial 
grid  are  calculated  in  subroutine  RAP. 

To  determine  whether  axial  rescaling  is  necessary,  densities  in  the 
cells  in  the  innermost  radial  shells  are  tested  m  BANGZ.  The  sum  of  the 
densities  in  the  first  six  cells,  SUMA,  and  the  sum  of  the  densities  in 
the  last  six  cells,  SUMB,  are  compared  with  the  sum  of  the  densities  in 
all  the  inner  cells,  SUMC.  If  one  of  the  ratios  SUMA/SUMC  or  SUMB/SUMC 
is  greater  than  0.00001  and  the  other  is  greater  than  0.000005,  rescaling 
of  the  axial  cells  by  a  factor  of  2  is  called  for.  If,  however,  one 
ratio  is  greater  than  0.00001  while  the  other  is  less  than  0.000005, 
indicating  that  the  ion  cloud  is  leaving  the  cell  space  at  one  end 
because  of  chemistry,  the  number  densities  are  shifted  by  one  cell  to 
reduce  the  problem.  If  rescaling  of  the  cell  space  has  been  called  for, 
due  to  diffusion  or  chemistry,  the  new  number  densities,  corresponding  to 
an  expanded  axial  grid,  are  calculated  in  subroutine  ZAP.  The  main 
program  suppresses  axial  rescaling  if  the  center  of  the  cell  space  of  the 
most  slowly  drifting  ion  has  drifted  a  distance  greater  than  the  length 


of  the  drift  tube. 
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VII.  INTEGRATION  OF  PREDICTOR-CORRECTOR  METHOD 
USING  CALI.  CAL2.  AND  S 


Calculation  of  number  densities  in  the  cell  space  of  each  ion  begins 
at  T  =  0  with  the  initial  values  mentioned  in  Section  IV.  For  convenience  we 
shall  describe  the  process  only  for  species  A.  The  appropriate  DT  is 
computed,  and  the  time  counter,  T,  is  incremented  by  this  factor.  The 
first  application  of  subroutine  S  in  each  iteration  computes  the  time 

J  /\  ■**  +  y-N  -C  4>U/\  -i  « -5  ♦  -I  n  1  i  m  f  A  nil  /»/%  11c*  1  «  4-  Vi  A  +  A/i 

ucuvaiivc  ui  luc;  liiiiiai  iiluu  uci  ucusj.  t)'  ,  rv ,  ivi  a.j.x  no  xn  cuv  tiu 

inner  radial  shells  and  stores  these  values  in  AP.  Then  CALI  uses  A  and 
AP  to  compute  the  predicted  value  of  the  number  density  at  time  T, 

denoted  by  PA ,  in  all  cells  in  the  three  inner  radial  shells  by  the 

first  order  Euler  method 


PA  =  A  +  AP  x  DT  .  (12) 

CALI  assumes  that  the  radial  distribution  of  any  ion  at  any  axial  position 
is  Gaussian  in  order  to  compute  P.A  for  cells  in  the  outer  radial  shell. 

The  technique  uses  the  PA's  for  the  three  cells  directly  inside  that  outer 
cell  and  a  relation  linking  segments  of  a  Gaussian  distribution,  as 
follows : 

AP(cell  4)  =  AP2(cell  3)/AP(cell  1)  ,  (13) 


where  the  inner  cell  is  cell  1  and  the  outer  cell  is  cell  4.  The 
time  derivative  of  the  predicted  value  of  A  is  next  calculated  by  S  for 
all  cells  in  the  three  inner  radial  shells,  and  the  values  are  placed  in 


__i _ *i  _  j _ 


Art'-,  rinaiiy,  lal^  calculates  tne  value  of  A  at  the  time  i,  replacing 
the  old  value  of  A.  For  cells  in  the  three  inner  radial  shells  the  second 
order  Euler  method  is  used. 


A  =  A  +  (1/2) * (APP  +  AP) *DT  , 


(14) 


where  we  have  written  a  FORTRAN  statement  defining  A  over  the  time 
interval  DT.  For  cells  in  the  outer  radial  shell,  the  new  A  is  also 
calculated  assuming  Gaussian  radial  distributions  by  the  formula 

A(cell  4)  =  A2(celi  3)/A(cell  1)  .  (15) 


1  n  i  T 


: _ n _ ^  i 


_ j  r* _  _ 


j  n  *  _ 


ax  xms  time  bainuk  ana  bain vjl  test  ror  tne  neea  lor  rescaling,  ana  kat  ana 
ZAP  are  used  for  this  purpose  as  needed.  A  new  DT  is  calculated,  T  is 
advanced,  and  calculations  begin  again  at  S. 

-  ~  9  --------  W 

21 


VIII.  CURRENT  CALCULATION  -  CURNT 


In  subroutine  CURNT,  the  currents  reaching  each  gate  are  calculated 
and  stored  in  AGAM,  BGAM,  and  CGAM,  etc.  The  routine  determines  the 
position  of  the  gate  in  each  cell  space  and  calculates  the  current  from 
the  densities  of  the  two  adjoining  cells  which  straddle  the  gate  at  the 
time  T.  That  is,  the  value  is  calculated  according  to  the  formula 


r 

z 


~  An 

-D  -r—  +  v  x  n 
z  Az  z 


(16) 


where  I\,  D_,  An/Az,  and  v_  are  the  values  at  the  gate  of  the  current, 
diffusion  coefficient,  space  rate  of  change  of  number  density,  and  drift 
velocity,  respectively,  in  the  axial  direction,  and  n  is  the  number  density 
at  the  gate.  The  value  of  n  in  Equation  (16)  is  determined  by  linear 
interpolation.  If  all  of  the  cells  in  the  cell  space  have  passed  the  gate, 
the  calculation  of  the  current  is  terminated  for  that  ion  at  that  parti¬ 
cular  gate.  When  the  cell  spaces  for  ions  A,  B,  and  C  have  all  passed  the 
final  gate,  the  entire  calculation  is  stopped  and  the  program  prints  output 
values . 


IX.  EFFECT  OF  CHANGING  MAXZ 

A  comparison  of  results  obtained  from  running  the  code  for  differing 
values  of  MAXZ  is  shown  in  Figure  2.  The  initial  parameters  of  the 
calculation  are  shown  on  the  figure.  The  use  of  MAXZ  =  63  or  greater  is 
seen  to  produce  essentially  identical  results,  except  that  the  use  of 
larger  values  of  MAXZ  yields  more  data  points  on  the  arrival  time  curve. 
The  run  time  for  the  code  also  increases  as  MAXZ  increases,  however. 

X.  COMPARISON  WITH  ANALYTICAL  CODE 
1* 

Woo  and  Wneaiton  have  developed  a  three-dimensional  analytical 
model  of  a  drift  tube  which  considers  drift,  diffusion,  and  one  reaction. 
Figure  3  shows  a  comparison  of  a  calculation  by  their  code  to  a  calcula¬ 
tion  made  using  MULTIGATE,  with  the  initial  conditions  for  each  calcula¬ 
tion  set  as  nearly  identical  as  possible.  The  arrival  curves  were 
normalized  to  produce  the  same  height  peaks,  but  are  otherwise  unchanged. 

The  agreement  could  be  improved  with  a  few-microsecond  shift  of  the  time 

_ 

References  may  be  found  on  page  26. ^ 
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axis,  but  the  curves  are  seen  to  be  in  agreement,  in  spite  of  the 
drastically  different  methods  used  to  calculate  them.  It  is  not  known 
if  the  shift  disagreement  is  due  to  assumptions  used  in  the  analytical 
model  or  to  truncation  effects  in  the  numerical  model. 

XI .  SUBROUTINE  PLOT 

Subroutine  PLOT  was  written  to  provide  a  rough  plot  of  the  three 
theoretical  current  distributions  on  a  single  computer  output  page.  The 
height  of  each  distribution  is  independently  scaled  with  ordinate  values 
referring  to  ion  A  only. 

PLOT  uses  two  other  subroutines,  SETUP  and  SORT.^  SETUP  is  used  in 
the  Y-scaling  and  returns  the  minimum  Y-value,  the  scale  factor,  and 
the  width,  defined  as  the  maximum  minus  the  minimum  value.  After  scaling, 
points  are  plotted  on  a  50  *  120  grid  using  the  letters,  "A,"  "B,"  and 
"C."  In  the  event  that  two  points  fall  in  the  same  position,  a  "D"  is 
printed. 


XII.  SUMMARY 

A  FORTRAN  IV  program  has  been  written  which  models  ion  transport 
and  chemistry  in  a  drift  tube.  While  the  program  presented  in  this 
report  is  concerned  with  three  ions,  the  program  could  be  extended  to  four 
or  more  ions,  if  desired.  Similarly,  the  program  could  accommodate  any 
conceivable  ion  distribution  at  the  beginning  of  the  ion  drift.  The 
program  can  therefore  be  used  to  model  the  ion  transport  and  chemistry 
of  any  conceivable  drift  tube  experiment.  It  can  also  be  used  to  check 
analytical  drift  tube  models  as  well  as  to  determine  the  conditions  under 
which  one-dimensional  drift  tube  models,  whether  analytical  or  numerical, 
can  be  safely  used. 
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APPENDIX 


LISTING  OF  PROGRAM  WITH  OUTPUT 

A  complete  listing  of  the  program  and  its  output  is  reproduced  on 
the  following  pages.  This  calculation  took  3.4  min  on  the  BRL  computer; 
it  could  be  expected  to  take  somewhat  less  than  1  min  on  a  CDC  6600 
machine. 
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*  SP 151  GEORGE  KELLER  394  X3335  MULTI  GATE  11 

S  MAXOI 5000 )L INES 

01  MENS  I  ON  A(  412  ),  API[  412  )  ,  PA(  412  )  ,TA(  1  03)  ,  AG  AM  ([412)  ,  APP  (412  ) 

DIMENSION  B(  412  ) ,  BPI[  412)  ,  PB(  412  > ,  TB(  1  03)  ,BGAH(  412  )  ,BPP  (412  ) 

DIMENSION  C(  412),  CP  1(412),  PC  (412  ),TC(103)  .CGAMI412)  ,CPP(412) 

DIMENSION  AGAMK  103  l,AGAM2(  103)  ,  AGAM3  (  103  )  ,  AGAM4  (  103  ) 

DIMENSION  BGAM1(103  ),BGAM2(  103)  ,BGAM3(103)  ,BGAM4  11103) 

DIMENSION  CGAM1( 103),CGAM2( 103) ,:GAM3(103) .CGAM4U03) 

DIMENSION  T A II 103 ) , I A2 d 03 ) , TA3 ( 1 03 ) , TA4 ( 1 03 ) 

DIMENSION  T B If 103>,TB2I 103 ) ,TB3 ( 103) ,TB4(103I 
DIMENSION  TCI! 103),TC2( 103 ) , TC3 ( 103 ) , TC4( 103 ) 

COMMON  MAXZ:,DR,DZ,Dir,PI 
P I  *3*  14 1592652)6 

C  REDUCED  MOBILITY  OF  SPECIES  A  IN  UNITS  OF  CM**2/V*SEC 
XK0A=3. 67 

C  REDUCED  MOBILITY  OF  SPECIES  B  IN  UNITS  OF  CM**2/V*SEC 
XKOB-1.85 

C  REDUCED  MOBILITY  OF  SPECIES  C  IN  UNITS  OF  CM**2/V*SEC 
XK0C*1. 65 

C  TOTAL  DRIFT  DISTANCE  IN  UNITS  OF  CM 
XL-7.5 

C  RADIUS  OF  INITIAL  ION  BURST  IN  RADIAL  DIRECTION  IN  UNITS  OF  ClM 
R0-.25 

C  LENGTH  OF  INITIAL  ION  BURST  IN  LONGITUDINAL  DIRECTION  IN  UNITS  OF  CM 
Z0*.3 

C  INITIAL  NUMBER  DENSITY  OF  SPECIES  A,  IN  ARBITRARY  UNITS 
DNA*1#  OE 12 

£C  INITIAL  NUMfliER  DENSITY  OF  SPECIES  B,  IN  ARBITRARY  UNITS 
DNBS 1. 0 

C  INITIAL  NUMBER  DENSITY  OF  SPECIES  C,  IN  ARBITRARY  UNITS 
DNC*l. 0 

C  GAS  TEMPERATURE,  IN  UNITS  OF  DEGREES  K 
TEMP-334. 

C  RATIO  OF  DRIFT  FIELD  TO  GAS  NUMBER  DENSITY,  IN  UNITS  OF  V*CM*CM 
EON-24. E-17 

C  AS  CURRENTLY  SET  UP,  THIS  PROGRAM  ASSUMES  THAT  THERE  ARE  TWO  GASES  IN  THE 

C  DRIFT  TUBE.  GAS  TYPE  1  IS  ASSUMED  TO  REACT  WITH  THE  IONS,  GAS  TYPE  2  IS  THE 

C  BUFFER  GAS.  CHANGES  TO  SUIT  OTHER  CONDITIONS  ARE  STRAIGHTFORWARD. 

C  PRESS JRE  OF  THE  REACTANT  GAS,  IN  UNITS  OF  TORR 
P 1  *  2  •  0 

C  PRESSURE  OF  THE  BUFFER  GAS,  IN  UNITS  OF  TORR 
P  2  *  0 . 0 


C 

RATE  CONSTANT  FOR 
YK AB  =  2  *  8E-30 

REACTION 

OF 

ION 

A 

TO 

FORM 

I  ONI 

B 

,  C,M**6/SEC 

c 

RATE  CONSTANT  FOR 
YKIBA-2. OE-1 4 

REACTION 

OF 

ION 

B 

TO 

FORM 

I  ONI 

A 

,  CM**3/SEC 

C 

RATE  CONSTANT  FOR 
YKBC=i.2E-29 

REACTION 

OF 

ION 

B 

TO 

FORM 

ION 

C 

,  CM* *6 /SEC 

C 

RATE  CONSTANT  FUR 

REACTION 

OF 

ION 

C 

TO 

FORM 

I  ONI 

8 

,  CM** 3 /SEC 

YKCB= 1. 58E-13 

FAC=( 3. S4E 16*273.0) /TEMP 

XN 1 -P l *  F  AC 

XN2-P2*F AC 

FAICTM  PUP2  )*273.0/(  760.  0* TEMP  ) 

XKA-XKOA/FACT 
XKIB-XKOb/FACT 
XKC-XKOC/FACT 
XNXN- XN 1  * ( X  N 1  +  X  N2 ) 

C  RAB  IS  THE  FREQUENCY  FOR  REACTION  OF  ION  A  TO  FORM  ION  B 
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RAB=YKAB*XMXN 
RBA=  YKEIA*(  XN  1L  +  XN2  ) 

RBC=YK8C*XNXN 
RCB=YKCB*(  XN1L+XN2) 

R  A  C  =  0 , 0 
RCA- 0*  0 

E«E0N*(  XNUXN2I 
C  THE  VELOCITIES  OF  THE  IONS 
V A=E*XK A 
V8=E*XKB 
V  C  =  E  *X  K.C 

FAClxTE:HP/l*  JL6E4 
D AR*  XK A*F AC  I 
D8R  =  XKEI*FAC1 
DCR=  X  KC*F  AC  1 

C  THE  PROGRAM  ASSUMES  THAT  THE  R-  AND  Z-D I RECT I  ON  DIFFUSION  COEFFICIENTS  ARE 
C  EQJAL.  THE  Z-DKRECT I  ON  COEFFICIENT  CAM  BE;  EXPECTED  TO  INCREASE  AT  HIGH 
C  VALUES  OF  E; / N • 

DAZ=DAR 
DBZ  =  DBR: 

OCZ  =  OCR 
OZ  *Z  0/310*  0 
OR*R  0/4* 5 

C  PRINTING  OF  INPUT  VALUES  AND  COMPUTED  PARAMETERS* 


PRINIT 

2101, 

P  1  ir  P  2 

PRINIT 

202, 

YKAB  t YKBA 

PRINIT 

2122  * 

YKBC  , YKCB 

PRINIT 

2105  , 

ZO„RO,,XL 

PRINIT 

2!06* 

TEMP 

PRINIT 

2107  » 

DN A, DNB, DSC 

P  R  I  Nl  T 

2110, 

XKOA, XKOB, XKOC 

PRINIT 

2!  15, 

E,  EON 

PRINIT 

200, 

VA„VB ,VC 

PRINT 

208, 

DAR, DBR, OCR 

PRINIT 

209, 

DA/ , DBZ  *  DC Z 

PRINT 

211, 

RAB, RB A 

PRINT 

231, 

RBC, RCB 

PRINT 

214, 

XN1, XN2 

C  SETTING  OF  INTERMEDIATE  GATE  POSITIONS* 

XL1=.2*XL 

XL2=  a  4*'XL 

XL3a;*6*XL 

XL4  =  .8*XL 

HZ*( ZO/ 2.0/DZ >**2 

HR*  I  RO/’DR-O*  !>)*’►  2 

C  N  DETERMINES  THE  NUMBER  OF  CELLS  IN  THE  Z-D I REC T I  ON ( MAXZ  = 2N  +  3 ) • 
N  =  40> 

MAXZ  =N*N+3 
X  =  N  ♦  1 
I  I  1=1 

1=0 

37  I F ( I.EUol I  GO  TO  10 

C  LIMITS  DT  TO  MAX  ALLOwlED  BY  DIFFUSION  OF  SPEC11ES  A 
DTAD= l« 0/ ( DAZ/DZ**2*DAR/DR**2 ) / 8* 0*3*0 
C  LIMITS  DT  To  MAX  ALLOWED  BY  DIFFUSION  OF  SPECIES  B 
D  T  BD=  1 3  0/  (  DH2/D2**2*DBR/DR**2)/R»0*3ia0 
C  LIMITS  ClT  TJ  MAX  ALLOWED  BY  DIFFUSION  OF  SPECIES  C 
0TCD=1. 0/1  DCi'/DZ**2*DCR/IH**2I  /B.0*3I,0 
C  LIMITS  DT  TO  MAX  ALLOWED  BY  THE  DRIFT  VELOCITIES  0 IF  THE  IONS 
DTY=  lo5*DZ/(  AMAXK  VA,  VB,  VC  )  +1*0  I 
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C  LIMITS  DT  TO  MAX  ALLOWED  BY  REACTIONS  AMONG  THE  SPECIES 
0TR*0.8*<  WO/AMAXK  RAB,rRBA,RACtRCA,RBCtRCB>  ) 

DT*  AMIN1(  DTADt  DTEJOt  DTCDtDTVtDTR  ) 

IF(DT.EQ.DTR)  I  *  1L 
10  GO  TO  ( 3 8 »  4 0 » 4 1 ) v 1 1 1 
C  SWITCHES  SET  TO  0  INITIALLY 
38  MAX A*0 
MAX  A1  =  0 
MAX A2=0 
MAX A3  *0 
MAX A4=0 
MAXB*0 
MAXB1 *0 
MAX  B2  *0 
MAXB3=0 
MAX  B4*0 
MAXC*0 
MAXC1 *0 
MAXC2*0 
MAXC3*0 
MAXC4*0 
L*0 
MB*0 
M1*0 
M2  *0 
M3*0 
M4*0 
K  A*0 
K 1  *  0 
K2*0 
K3*0 
K4*0 
LC*0 
LI  *0 
L2*0 
L3*0 
L4*0 
I  A,A*0 
I  A A  1*0 
IAA2»0 
I  AA3=:0 
I  A A 4  =  0 
IBB*  Cl 
I BB 1*0 
I  BiB2~0 
I BB3-0 
IBB 4*0 
ICC*0 
ICC1-0 
1CC2-0 
ICC 3*0 
ICC4-0 
I  NIGA«0 
INGA 1*0 
I NIGA2  =  0 
I  NIG  A  3*0 
I NG A4*0 
I NGB  =  0 
I  M  G  B 1L  -  0 
I N  G  B  2  =  0 
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I N  G  B  3  =  0 
I NGB4  =  0 
I  NGC:=  0 
ingc:l=o 

I NGC2=0 
I NGC3=0 
I NGC4-0 

C  SUMS  UNIT  I  ALLY  SET  TO  0 
AGAM ( MAXZ ) =0* 0 
AGAM1(MAXZ)*0.0 
AGAM2 ( MAXZ )=Q.O 
AGAM 3 ( MAXZ )=0.0 
AGAM4 ( MAXZ 1=0.0 
BGAM( MAXZ >=0.0 
BGAMK  MAXZ  1*0.0 
BGAM2IMAXZ )*0.0 
BGAM3( MAXZ ) *0.0 
BGAM4I MAXZ) *0.0 
CGAM( MAXZ )*0.0 
CGAMK  MAXZ  )  =  0.0 
CGAM2 ( MAX  Z )  *0.  0 
CGAM3 ( MAXZ )* 0.0 
CGAM4IMAXZ )  =  0. 0 

C  CAL.CJL  AT  I  ON  OF  INITIAL  OEMSITY  DISTRIBUTIONS. 

DO  6  K*  I  *  4 
XK*  K 

DO  6  J=  1 1  M AXZ 
XJ*J 
L  :=  L  ♦  1 

A ( L ) =EXP( -( XJ-X )**2/HZ-< XK-0.5 )**2/HR) 

I  IF  (  J  •  EQ  >  1  )  A  (  L  )  =  0  >  0 
C(L  )  =  A(  L  )  ♦DNC/DNI A 
6  B(L)  =A(  L)*DNB/DNlA 

C  DAff  DB • DC  REPRESENT  TOTAL  DISTANCE  TRAVELED  BY  A,B,C. 

DA=0,0 
DIB*0o  0 
DC=0. 0 
T  *  0  •  0 

DO  100  JKL*  L  « 1000 
T  =  T  ♦  D  T 
DA*DA*V  A*DT 
OB=DB*VB*DT 
DC*DC+VC*DT 
L  *0 

C  AP  IS  THE  DERIVATIVE  AT  THE  BEGINNING  OF  THE  INTERVAL. 

CALL  S  (  A,  B>,  C  t  AP  ,  D  AZ  ,  D  AR  ,  R  AB  ,  RB  A  ,  R  AC  ,  RC  A  ,  DA  ,  Dfll ,  DC  ) 

CALL  S ( B*  C, A, BP, DBZ ,DBR, RBC.RCB ,RBA,RAB,DB,DC ,DA ) 

CALL  S(C,A,B,CP,  L)CZ  ,  OCR,  RC A , R A C  ,  R C B . RBC * DC  t DA  * DB I 
C  C  CALI-FIRST  - OiRDER  PREDICTOR  SCHEME  (  P  A  ,  PB  ,  PC -PRECI 1  C  TED  VALUES). 
CALL  CALK  A,  AP,  PA) 

CALL  CALK  B,  BP.PB) 

CALL  CAL  1 ( C*  CP*  PC ) 

CALL  S  I  P  A ,  P  B  *  P  C  ,  A  P  P  ,  D  A  Z  ,  0  A  R  ,  R  A  B  ,  R  B  A  ,  R  A  C  ,  R  C  A  ,  D  A  .  D  B  #  DC  ) 

CALL  S( PB, PC,PA,BPP,DBZ,DBR,RBC , RC B . RB A , RAB ,0B , DC , DA ) 

CALL  S(PC*PA,PB,CPPtDCZ,DCR,RCA,RAC,RCB,RBC,DC,DA,DB) 

C  CAIL2  SECOND  ORDER  CORRECTOR 
CALL  CAL2I A, APP, AP) 

CALL  CAL  2 1 B, DPP  *  BP ) 

CALL  CAL2IC.CPP*CP) 

C  TEST  FOR  RADIAL  RESCALING* 
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CALL  BANGRII  AyMAXZySING) 

CALL  BANGRII  By MAXZ,  TING) 

CALL  BANGRII  Cy  MAXZ yZ  ING) 

IFlSING*TING*ZING.EQ.O.O)  GOTO  40 
C  Ay  By  C  RESCALED  IN  R-DIRECTION  BY  RESETTING  DR  TO  l.l*DR. 


CALL  RAP I  Ay  HAXZ ) 

CALL  RAPI  B«  MAXZ  ) 

CALL  RAP! Cy MAXZ) 

DR*DR*l«l 

111*2! 

GOTO  37 

40  CONTINUE 

C  TEST  FOR  Z  RESCALING 

CALL  BANGZIl  Ay  HAXZ»SING«DA«DZ»KA  ) 

CALL  BANGZIlByMAXZyTINGyDByDZyMS  ) 

CALL  BANGZIl  C.MAXZyZ  INGyDCyOZyLC  ) 

IFIS1NG*TI'IG+ZING.E0.0.0)  GOTO  41 

C  IF  SLOWEST  ION  HAS  HAD  CHANCE  TO  PASS  OUT  OF  TUBE  ( NEGLECTING  DIFFUS IIONII  • 

C  STOP  Z  RESCALING 

IF(C!AHKNlllVAyVByVC))*T).GT.XL)  GO  TO  41 

C  A« By C  RESCALED  IN  Z-DIRECTION  BY  RESETTING  DZ  TO  2.0*DZ. 

C  IT  IS  ABSOLUTELY  ESSENTIAL  THAT  THE  FACTOR  BY  WHICH  Z  IS  RESCALED  BE  LEFT  AT  2 
CALL  ZAP!  Air  HAXZ) 

CALL  ZAP!  B«  HAXZ) 

CALL  Z  API  C»  HAXZ) 

DZ-DZ>DZ 
111*31 
GOTO  37 

41  CONTINUE 

C  CURRENTS 

CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 

C  TEST  TO  SEE  WHETHER  DISTRIBUTIONS  HAVE  PASSED  THE  LAST  GATE. 

C  CALCULATION  STOPS  WHEN  THE  3  DISTRIBUTIONS  HAVE  PASSED  THE  FINAL  GATE. 
IFIINIGA.E0.1.AND,  INGB.EO.l.AND.  INGC.EQ.I)  GOTO  101 

100  CONTINUE 

C  OUTPUT  VALUES  PRINTED  FOR  FINAL  AND'  4  INTERMEDIATE  GATES. 

101  CALL  OUTPI AGAH1 y  TA lf BGAMly TB l y CGAHl * TC 1 ) 

CALL  0UTPIAGAM2,  T A2 , BGAH2 y TB2 y CGAM2 y TC2 ) 

CALL  0UTPIAGAM3,  TA 3 • BGAM3 y TB 3 y CGAM3 y TC3 ) 


CALCULATED  AT  5  GATES  FOR  EACH  SPECIES. 

CURNT I  Ay  MAX  Ay  N  y  VAy  T  «XL  y  DAZ  y  TA  y  I AA  »KAy AGAMy INGA  yDA) 

CURNT  II  Ay  MAXAly  Ny  VAy  Ty  XL1  y  DAZ.TA1  y  I  AA1  yKl y AGAM1  y  INGA1  y  DA) 
CURNT I  Ay MAX A2y  Ny  VAy  Ty  XL2V  DAZ  y  TA2y I AA2  y K2  « AGAM2  y INGA2  yDA) 
CURNT II  AyHAXA3yN«VA»TyXL3yDAZy  TA3»  I AA3  VK3 «AGAH3  y  INGA3  y  DA) 
CURNT  II  Ay  MAX  A4V  Ny  VA«  Ty  XL4y  DAZ.  TA4y  I A  A4  y  K4  y  AGAH4  #  I NGA4  v  DA) 
CURNT  II  ByMAXByNy VByT  yXLyDBZ  yTly  IBByHByBGAM*  INGByDB) 

CURNT  IE  fiyHAX  Bly  Ny  VBy  Ty  XL1  yDBZt  TBly  IBB1  y  Ml  yBGAHl  v  1NGB1 1 DB) 
C  JRNTII  By  HAXB2*  N«  VBv  Ty  XL2y  DBZ  t  TB2V IBB2  y  M2«BGAM2  y  ING82  *DB ) 
CURNT II  ByHAXiB3yNv  VBy  TyXL3yDBZ»TB3yIBB3yM3yBGAM3yINGB3yDB) 
CURNT II  ByHAXB4»  Ny  VBy  Ty  XL 4 y  DBZ  •  TB4,  IBB4» M4f  BGAH4  •  I NGB4  y  OB) 
CURNT II C  yHAXCy  N  •  VCy  T  y  XL  y  DCZ  *TC  y  ICC  y  LCy  CGAMy  I NGC  y  DC  ) 

CURNT1I  CyHAXCly  Ny  VCy  TyXLlyDCZy  TCly  ICC1  yLl  tCGAMly  INGC1  yDC) 
CURNTII  Cy  HAXC2y  N,  VC  y  Ty  XL2  y  DCZ  y  TC2  y  ICC2  y  L2  9 CGAM2  v  INGC2  »  DC ) 
CURNTI C *  HAXC3V  Ny  VC  v  Ty  XL3 »DCZ  y  TC3 y ICC3  y  L3 y CGAM3 y I NGC3  y  DC ) 
CURNTI[CyHAXC4»NyVCyTyXL4yOCZyTC4,  ICC4  ,L4f  CGAH4 ,  INGC4 y  DC) 


CALL  OUTPI AGAM4,  TA4,BGAM4, TB4, CGAM4 , TC4 ) 

CALL  OUTPI AGAMy TAyBGAMyTByCGAMy TC ) 

STOP 

200  FORMAT!  1.H  ,  30X,  3HVA  *,  1PE10.  2  y  10X,  3HVB*  y  1  PE  10.2  y  1  OX  ,3HVC*  y  1  PE  10.  2  ) 

201  FORMAT!  IHly  30X, 3HP 1 - y F 7. 3 • 1 3 X y 3 HP 2 * y F 6 • 4 ) 

202  FORMAT!  IH  y  28X  y  5HYK  AB  =  ,  IPE  10.  2  ,  8X  y  5HYKBA**  JLPE3L  0.  2  ) 

205  FORMAT!  IH  *  30X,  3HZ0  =  ,F6.4,  14Xy  3  HR  0  *tF  6.  4, 1 4  X,3HXL*y  IF  9.  5) 
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206  FORMAT C 1H  , 28X, 5HTEMP-, F5 . 1 1 

207  FORMATt  IH  ,  29X,  4HDN  A*,  1PE  10.  2*  9X,  4H0NB- , 1PIE10.  2  ,  9X  ,4H0NC»  » 1  PE10.  2  !► 

208  FORMATt  !LH  ,  29X,  4H0AR =  ,  1PE  10,  2, 9X,  4HDUR*  ,  1PIE10.  2 , 9X  * 4H0CR*  1 1PE10. 2  I 

209  FORMATt  LH  ,29X,4HDAZ»,  1PE  10.  2 1 9X, 4HDBZ*  , 1PIE 10, 2  * 9X , 4H0CZ-  ,  IPE1C1.  2  II 

210  FORMATt  LH  ,  28X,  5HXX  OA*,  F7  •  5 ,  11 X  ,  5HXK0B*  ,  F7<,  5 , 1 1  X  ,5HXK0C*  t  F7.  5  1 

211  FORMATt  LH  ,  29X,  4HRAB*,  1PE  10*  2,  9X,4HRBA*  ,  1PIE10.  2  1 
222  FO  RMATt  LH  ,  28X,  5HYKBC»,  IP  E 10. 2 „  8X ,  5HYKCIJ* •  1PE10.  2  1 
231  FO  RMATt  LH  ,  29X,  4HRBC*t  1PE  10.  2. 9X,  4HRCB*,  IPIE10.2  1 
216  FORMATtLH  ,  29X,  4HXN  1  *,  1PE  10.  2, 9X,  4HXN2*  , 1PIE  10.  2  1 
215  FORMATtLH  ,  3  IX,  2HE*,  1PE 1 0. 2 1  9X „  4HE0N *,  1IPE1  0. 2) 

ENID 

C  BANGR 

SUIBROJT  I  ME  BANGR  t  A,  MAXZ  ,  S  ING  » 

C  BANGR  DETERMINES  WHETHER  RESCALING  IS  NECESSARY  IN  THE  R-OI RECT I0N«. 
DIMENSION  At  1 ) 

C  L-DES I GNATES  INNER  RADIAL  ZONE. 

L»MAXZ/2 

C  K-DESIGNATES  OUTER  RADIAL  ZONE. 

K-L+MAXZ+MAXZ+MAXZ 

SING»0.(> 

C  RADIAL  RESCALING  CRITERIA. 

IF  t ACKl/AtL l.GT.0.90)  SING- 1.0 

RETURN 

END 

C-  IIANGZ 

SUIBROJT  INE  BANGZt  A,  MAXZ,  S  ING,DA,DZ,  I  API 
C  BANGZ  DETERMINES  WHETHER  RESCALING  IS  NECESSARY  IN  THE  Z-DIRECTION, 

C  IF  THE  NEED  OF  Z  RESCALING  IS  FOUND  AT  ONLY  ONIE  END  OF  THE  ION  CLOUD,  THE 
C  ENTIRE  CLOUD  IS  SHIFTED  ONE  ZONE  IN  THE  l  DIRECTION. 

DIMENSION  At  11 

C  SUMA-SJM  OF  DENSITIES  IN  63  CELLS. 

SUIMA«0.0 
DO  10  J-1,6 

10  SUMA«S JMA+At J ) 

S'JIMB»0.0 
MNI6-MAX2-6 
DO  20  J*MM6 , MAXZ 
20  SUIMB»SUMB*Al  Jl 
MM7*MM6  -1 
SUIMCss0.  0 
DO  30  Js 7, MM7 
30  SUIMCasSUMC*At  J  1 

SUIMC«SUMC+SUMB+SUMA 

RAT  A«SUMA/SUMC 

RATB«SUMB/SUMC 

ISL-0 

IS,R«0 

BZ-0,,00001 
BZZ  *0. 000005 
IF tRATA.GT.BZ 1  ISL*1 
IF tRATB.GT.BZ)  ISR*ISL*1 
S I NG« 1. 0 

IF  t  ISR.IE0.2)  RETURN 
IF  t  ISL.IEO.O)  GOTO  60 
MM*M AXZ — 1 

IF  (  RATB,,GT.  BZZ  I  RETURN 
DO  50  K  =  1 , 4 
MHM*K*MAXZ 
DO  50  J 3  L , MM 
L*  MMM-J 
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50  ACL  +  il-ACL) 

IAP*IAPU 
DA»DA-D  l 
SI  MIG*  0.0 
RET  URN 

60  IFt  ISR.ECI.0)  GOTO  90 
IFC RATA. GT.B ZZ  >  RETURN 
MN-MAKZ-l 
00  71  K«l,4 

MMM*  M  AX  Z  **  C  X  - 1 1 
00  70  J«2,MM 
L«JI*MMM 

70  ACL.  )*A(  L^l  | 

MZ*‘K*NAXZ 

71  ACMZI* A(MZ-l)/2,0 
IAP'»IAP-lL 
DA«DA*DZ 

90  SIMG»0.0 
REIURM 
END 

C  CALI 

SUDROUT IME  CAL  1C A,B,C) 

C  CALI  IS  A  FIRST  ORDER  PREDICTOR  SCHEME. 

C  A  DENOTES  THE  INITIAL  VALUE,  B  THIE  DERIVATIVE,  AND  C  THE  VALUE  PREDICT 
DIMENSION  ACDiBI  l),CCl> 

COMMON  MAXZ , DR, DZ ,DT ,PI 
L  *  Cl 

DO  10  K  =  l»  3 
DO  10  JS1L , MAXZ 
L»L.*1 

10  CCL)»ACL1I*B(L)  *DT 
DO  30  J*2, MAXZ 
Ll-'J 

L3*‘MAi;Z+MAXZ*J 
L4»>  L 3+MAXZ 
CCL4»«0.0 

30  IFC  CCLl  l,(GT.  l.OE-lOO)  CC  LA)  *C  C  (  L3  )  /C  ( L 1  \)*C  ( L3 ) 

RETURN 

END 

C  CAL2 

SUBROUTINE  CAL2(A,B,C) 

C  CAL 2  IS;  A  SECOND  ORDER  CORRECTOR 
DIMENSION  At  II, BC  1),C(1) 

COMMON  M  A  X  Z  ,  DR  ,  DZ  ,  D  T  ,  P  I 

DTT  =  D1  V2..0 

L*CI 

DO  10  K  - 1 , 3 
DO  10  JSL,MAXZ 
L  =  L>1 

10  AC  L  )*  AC  L  H  +C  B(L I ♦C € L I I*DTT 
DO  30  J=2,MAXZ 
L  1*>  J 

L3-  MAXZ+MAXZ  *J 
L4*  L 21  +  MAXZ 
AC  1.4)  «030 

30  IFC  ACL1UGT.  l.QE-100)  AC  L4  )  =C  A  (  L3  I /A  C  1. 1 ))  ♦  A  C  L3 ) 

RETURN 

END 


C 


SUBROUTINE  CURNTC  A, M AX, N , V A , T, XL , DAZ , T A , I  A A  ,  I AP , AGAM , I NG, DA ) 


CURNT 
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C  CURNT  GENERATES  THE  CURRENTS  AT  EACH  GATE, 

DIMENSION  At  ll.TAt 1 ), AGAMt  1) 

COMMON  MAXZ , DR , DZ ,DT,  PI 
IF( VA.EQ.0,01  RETURN 

C  I NlG*  1  INDICATES  THE  DISTRIBUTION  HAS  PASSED  THE  GATE, 

C  IAA  IS  THE  COUNTER  FOR  STORING  THIE  CURRENTS  AND  CORRESPONDING  TIMES  IN 
C  THE  AGAM  AND  TA  ARRAYS  RESPECT  II  VEIL  Y. 

IF<  ING«EQ,.1.0R,  IAA.GE.NAXZ-I)  GOTO  60 
C  XA  IS  THE  DISTANCE  FROM  THE  LEFTMOST  CELL  OF  THE  CELL  SPACE  TO  THE  GRI 
X  A*  «  XL4 FLOAT  I  N*  1I*DZ~DA  )  /DZ 
I  A=XA 
I  NGS!0 

I Ff  IA.GE.MAXZ  )  GOTO  50 
I F(  MAX, EQ» 0 )  I A  A  *  0 
IF!  MAX,,  EQ,  01  IAP*0 
Ml  AX*  I 

C  TEST  TO  DETERMINE  WHETHER  DISTRIBUTION  HAS  PASSED  THE  GATE, 

IF< I A,L  T, l ) G0T060 
L.2*  I A+MAX2 
AL2*A( L 2  I 
A J  *  A ( 1A) 

I F( AL2«  GT ,0, 0, AND, AJ.  GT •  AL 2 )  GOTO  AO 
GOTO  70 

AO  T  HE  atX  A- FLO  AT  (  I A  9 
0THE*1, O-THE 
I  X  *  I A 
A I  A* AJ 

A  I A 1 s  A (  I A*1 ) 

MOLD=OTHE*A(  IX  I +THE*A( !«♦! ) 

HOLD* -0 AZ  * C  A  I A1 - A 1  A ) /DZ ♦ VA*HOLD 
AGAM(  MAXZ  )*AGAM(  MAXZ  l«-HOLD*DT 
£  I  FI  I  A , E  Q , IAP  I GO  T 070 

I AP» I  A 
IAA»IAA+l 
AGAMt  I  A A )  =HQLD 
TA(  1[AA»=T 
50  RETURN 
60  I NG»1 
70  RETURN 
END 

C  INTER 

FUNCTION  INTER! X) 

C  INTER  INSURES  ROUNDING  DOWN  IN  ALL  CASES 
INTER*! X41000.0 ) 

INTER* INTER- 1000 

RETURN 

END 

C  OUT  P 

SUBROUT INE  OUTP < AGAM, TA,  BGAM, TB  ,  CGAM*  TC ) 

C  OUTP  PRINTS  THE  OUTPUT  VALUES.  OUTPUT  CONSISTS  OF  CURRENTS  AND  CORRES 
C  PONOING  TIMES  FOR  3  IONS  AT  5  GATES,  AS  WELL  AS  MEAN  DRIFT  TIMES,  SUMS 
C  OF  CURRENTS,  AND  RATIOS  OF  SUMS. 

DIMENSION  AGAM<  1 1  , TA ( 1 1 , BGAM( 1)  ,TB(l) , CGAMt II , T  C 1 1 1 

COMMON  MAXZ , DR, OZ, DT,  PI 

COMMON  /BAR/  TI,T2,T3 

PRINT  1L0 

DO  3  J -I, MAXZ 

IF(1A( J  )4TB( JI4TCC JI.EQ.0,0)  GOTO  4 
3  PRINT  1,  AGAMt  J  )  ,TA(  Jl)  ,BGAM(J  ),  TB  (  J  )  ,CGAM(J>  ,TC  (  J) 

PRINT  5,  AGAMt  MAXZ  ) 
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PRINT  7,  BGAM(MAXZ) 

PRINT  9,  CGAMIMAXZ) 

T  1*TBAR(  AGANt  TAt  SUN  II 
T2*TBAR(  BGAN*TBf  SUN  !> 

T3*TBAR( CGANf  TC«  SUM  I 

PRINT  666*  T I*  T2»  T3 

CALL  PL0T<TAf AGAM,TB»BGAMf TCtCGAM) 

RETURN 

5  FORMAT ( 1H  » /, 20X, 5HSUMA*. IPE12. 4) 

7  FORMAT <  1H  ,  20Xt  5HSUMB*.  1PIE  12.4 ) 

9  FORMAT  (  1 H  f  20X.  5HSUMC*  v  1PIE  12.  4 ) 

1  FORMAT (  IM  t  1P6E20. 4) 

10  FORMAT (  1 H 1  *  1 4X «  lHAf  117X. 4HTIME  •  1  8X.  1HB  •  17Xf4HTI  ME  »18Xf  1HC,  ]L7Xt 
14HTIME) 

665  FORMAT  (  1H  ,20X,  3HTA*f  1PE  12. 4. 5X  ,3HTB*  1 1 PE12. 4. 5X „ 3HTC* .  1PIE12. 4 ) 

ENID 

C  PLOT 

SUBROUT INE  PLOTIXlfVlt K2t Y2tX3f Y3) 

C  PLOT  PILOTS  3  ION  ARRIVAL  SPECTRA  VERSUS  TIME.  IT  IS  SET  UP  FOR  132  COL.  PRINT 
C  MAGIC  NUMBER  186  REFERS  TO  MAXIMUM  NUMBERS  OF  CURRENT/TIME  PAIRS 
C  MAGIC  NUMBER  62  REFERS  TO  MAXIMUM  NUMBER  OF  PAIRS  PER  ION  SPECIES 
INTEGER  BLANKt  CAt CB«  CC«  CD 

DIMENSION  Xim,X2<  UiXBUItYK  l)»Y2li)  fV3Cll 
DIMENSION  LINE!  I20)f M<  1L86),N(2*  186)  *RMAN«6)  ,R0RDtl3) 

DATA  CAfCBtCCtCDf  BLANK/ 1 HA « 1HB  * 1HC  » 1HD» 1H  / 

OATA  LINE/120*1H  / 

DO  2  I* 1 f 186 
MI I)«0 

2  Nllfl)*0 

CALL  SETUPlTlt YISf Y1SC.Y1D) 

CALL  SETUP! Y2fY2S»Y2SC,Y2D) 

CALL  SETUP! Y3,Y3SfY3SCfY3D) 

A*«1.E50 

B*1.E50 

L«0 

DO  6  I*  If  62 
L*IL  +  1 

N!lfL)*!Yl!  I  )-Y  IS  )*YlSC+0«  5 
N!2tL)sCA 

IFIINI  i.L  ).EQ.O)  GO  TO  4 
XX*X1 ( I  ) 

IFIIXX.GT.A)  A=X  X 
IFIIXX.LT.B)  B-XX 

4  L  *  L  ♦  l 

N! l»L)-(Y2t I )-Y2S)*Y2SC>0.5 
N  (  2 1  L  ) s  C  IB 

IF!  N(  I*  L  )  .  E 0 * 0  )  GO  TO  5 
X  X  *  X  2  {  I  ) 

IFIIXX.GT.A)  A=  XX 
IFIIXX.LT.B)  8=XX 

5  L*L+1 

N !  1L  t  L  )  =  (  Y 3 !  I  ) -Y  3S  )  *Y  3SC  ♦0.  5 
N! 2 «  L ) »CC 

IFIINI  I,  L  ).EQ.O)  GO  TO  6 
X  X  «  X  3  (  I  ) 

IFIIXX.GT.A)  A*X  X 
IFIIXX.LT.B)  B=XX 

6  CONTINUE 
X  1 D  =  A  -  B 
XlS=B-XlD/20. 
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XID«XID*1.1 

C  REPLACE  120  WITH  60  FOR  72  COLUMN  OUTPUT 
X  ISC* 120* /X 1 0 
L*0 

DO  7  1*1*62 
L  »  L  ♦  l 

M!IU*(Xim-XlS  I  +  X1SC  +  0.5 
L  =  L  ♦  1 

M II  L  )  *  (  X  2  (  I )-Xl$)*XlSC+0.5 
L  *  L  ♦  1 

MI|L»*(X3(I)-X1SI*X1SC^0,5 

7  CONTINUE 

DO  8  1=1*186 

I  IF  (  N (  1  *  IUGT.51I  N(  l*  I  I  *N(  1  *  I  I  "  50 

8  NU,n*51-NI  1,1  I 
CALL  S0RT<M,NI 
DO  9  1*1*13 

9  RORDI  I  >  =X  lS-fFLOAT (  I  - 1 II *X  1LD/ 12. 

DO  10  1=1*6 

10  RMAN(  I  )  =Y  IS  ♦FLOAT!  6>-  I  II  *Y  10/5* 

1  =  0 
LIL*  1 

111  I  *  I  ♦  l 

I  IF  (  N  (  1*  D.NE.O)  GO  TO  12 
GO  TO  11 

12  L*0 
PRINT  101 

13  L  *  L  ♦  1 
K*L-1 

K!X*(K/10)*10 
1  IF  (  L*  E0«  5  1 )  GO  TO  20 
£  21  I  IF  (  M  (  II.NE.O)  GO  TO  22 

I  IF  (  I  •  GT  o  1 8 6  )  GO  TO  23 
1*1+1 
GO  TO  21 

22  IF ( L • E  Q  >  N ( 1,1))  GO  TO  15 

23  IFIK.EQ.KX)  GO  TO  14 
PRINT  103 

GO  TO  13 

14  PRINT  104,  RMIAN(LL)! 

L  !L  =  L  L  ♦  1 

GO  TO  13 

15  J*M( l ) 

J  J*L I NE ( J  ) 

L INEI  J ) =N( 2,  l) 

IF(JJ*NE*8l  ANIK  )  LIMECJ)=CD 
1*1+1 

IF(L.EO»N(  1*  IM  GO  TO  15 
IF(K.EOoKK)  GO  TO  17 

C  REPLACE  120  WITH  60  FOR  72  COLUMN  OUTPUT 
PRINT  106,  (LINE!  J  II  •  J  =  1,120) 

GO  TO  18 

C  REPLACE  120  WITH  60  FOR  72  COLUMN  OUTPUT 
17  PRINT  105,RMAN<LL )* < L I NE (J ) , J* 1 , 120) 
LL=LL*i 

10  DO  19  J=l» 120 

19  LINE! J)=BLANK 
GO  TO  15 

20  PRINT  102,RMAN(LL)«(R0RDlJ),J=i,13i 
RETURN 
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101  FORMAT!  57H1THE0RETICAL  CURVES  -  THE  Y  AXIS  APPLIES  TO  ION  SPECIES 
LA) 

102  FORMAT (  E9.2,  12(  10H* - )  v  lH+»/»3X9i3EI0«2  I 

103  FORMAT  !  9X  »  1IH  I) 

104  FORMAT! 69. 2* 1H*  ) 

105  FORMAT! E9. 2, 1H*,120A1) 

106  FORMAT! 9X  ,  lH  I , 120A1 ) 

END 

c  RAP 

SUBROUTINE  IRAP!A,MAXZ) 

C  RAP  RESCALES  DENSITIES  IN  THE  R  DIRECTION. 

DIMENSION  At  1) 

DO  58  J*2,MAXZ 
L2* J+MAXZ 
L  3*L2^M AXZ 
L4*L3*MAXZ 
AJ-A(  J| 

AL3* A! L  3 ) 

I F( AL3. GT.O. 0. AND. AJ.GT. AL3 )  GOTO  56 

A! J  > *0. 0 

A!  L2 1*0,0 

A !  L  3  ) »  0  j  0 

A!L4)*0  0 

GOTO  58 

56  ALP*ALOG! AJ/AL3I/24.0 
DO  55  K  =  ;if  4 
L4*  J  +  l!  K-l  )*MAXZ 
HA* 1. 2 1 *FLOAT! X ♦X-l 1**2  -1.0 
&  55  A!L4)*AJ*EXIP(-ALP»HA) 

58  CONTINUE 
RETURN 
END 

C  S 

SUBROUTINE  S!  A,  B ,  C ,  AP»  D  AZ  ,  DAR  ,  RAB  »RBA  ,  RAC  .IRCA  ,DA,DB  .  DC  I 
C  S  RETURNS  THE  DERIVATIVE  IN  AP 

DIMENSION  A!  L),AP(l),&!  I  I  ,  C  ( 11  *  GAR  1 103 )  fGARlP  ( 1  03  ) 

COMMON  MAXZ, DR, DZ,DT,PI 

C  XLC  MEASURES  THE  RELATIVE  DISTANCE  OF  C  FROM  A  IN  UNITS  OF  OZ, 

X LC* !  DA-DC ) / DZ 
LPC* I  NT cR( XLC ) 

T HEC*XLC-FLOAT<  L PC ) 

C  X  LB  MEASURES  THE  RELATIVE  DISTANCE  OF  ES  FROM  A  IN  UNITS  OF  DZ. 

X  LB* ! DA-DB )  /DZ 
LPB* I  NT £R ( XL  8 ) 

THEB*XL H-FLO AT!  LPB) 

OTHEB*1jO-ThEB 
OTHEC* 1.0-THEC 
L  =  1 

LR*M AXZ 
DO  12  K  *  1L  ,  3 
X  K*K 

XKP=K  +  K-]L 

T  PKDZ  *2 , Q*P I *XK  *DZ 
ARA«P|[*XKP*DR**2 
VOL* ARA*DZ 
DAZDA  =  OAZ/DZ  *AR  A 
DART P* DAR *T PKDZ 
HAP*0. 0 
G AZ  =  0. 0 

DO  12  J  * 1 , M  A  X  Z 
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LM1  =  L 
L  *L  ♦  1 
LR=LR+1 

IF(K.GT.i)  GOTO  10 
GAR<  J )  =  0®  0 
10  GAZ  P  =  GAZ 
HA=0.0 

I  F ( J*LT »MAXZ )  HA  =  A( L ) 

GAZ=-OAZDA»<  HA-HAP) 

HAR  =  A ( L  R ) 

GARP ( J ) =GAR (  J ) 

GAR( J )*-DARTP*(  HAR-HAP ) 

IF(J.EQjl)  GOTO  14 
JP8* J+LPB 

KPB=JPB*MAXZ*(X-i  I 
JPCS J+LPC 

KPC*JPC  +  MAXZM X-l  ) 

R*RAB+R AC 
CHB=0.0 

IF( JPB.LT.HAXZ.AND. JPB.GT.O  )  CHB*RBA* (B( KPB ) *OTHEB 
1*B( KPB+l >  *THEB ) 

CHC- 0*  0 

IF( JPC*LT.MAXZ.AND. JPC.GT.O  )  CHC*RCA* ( C ( KPC I *OTHEC 
1*C<KPC+1)*THEC) 

API LM1 ) *  -(HAP*R-CHB-CHC«-(GAR ( J  )-GARP( J ) ♦GAZ-GAZP ) /VOU 
14  CONTINUE 
HAP=HA 
12  CONTINUE 
RETURN 
END 
C 

^  SUBROUTINE  SETUP ( X, B, Yi SC , Y ID ) 

°  C  SETUP  DETERMINES  SCALE  FACTORS  FOR  PLOTTING  DIFFERENT  ION  SPECIES 
C  MAGIC  NUMBER  62  REFERS  TO  MAXIMUM  NUMBER  OF  PAIRS  PER  ION  SPECIES 
DIMENSION  X( 1) 

A  =  - 1*  E50 
B= 1 • E50 
DO  1  Ls 1  *  62 
XX*X ( L I 

IFIXX.GT.A)  a*xx 
IFIXX.LT.B)  B=X  X 
1  CONTINUE 
Y 1 D=  A-B 
Y 1 SC  =  50„ 0/Y l D 
RETURN 
END 
C 

SUBROUTINE  SORT(M,N) 

C  SORT  INTERSORTS  CURRENT/TIME  PAIRS  PRIOR  TO  THEIR  BEING  PLOTTEO 
C  MAGIC  NUMBER  186  REFERS  TO  MAXIMUM  NUMBERS  OF  CURRENT/TIME  PAIRS 
DIMENSION  MI  186), Nl 2,186) 

MM= 1 8  6 

1  MM=  MM/2 

IFIMM.EO.O)  RETURN 
K= 1 86-MM 
J  =  1 

2  I=J 

3  IM*l+MM 

I  F  I  N(  1,  I  )*LE.N(  1,  IM) )  GO  TO  4 
I  T  =  M  (  I  ) 


SETUP 


SORT 
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Mil  I  1  =  M I  IM) 

Mil  I M  )  *  I  T 
DO  5  L ss  1  *  2 
IT«MlLt  I) 

Nil  L,  I  )  ==  N l  IL,  IM) 

5  NHL,  IM  11  =  IT 
l *  I —MM 

I  FI  I.GE  ,1  )  GO  TO  3 
4  J»J*1 

IFl J.GT > K I  GO  TO  I 

GO  TO  2 

END 

C  TBAR 

FUNCTION  TBAiRlXNfTI,  SUM  ) 

C  TBAR  CALCULATES  THE  MEAN  DRIFT  TIMEo 
DIMENSION  X  N  l  1 )  ,  T  1 1  L) 

COMMON  MAXZ, DR,DZ,DT,  PI 

Ml  =  M  A  X  Z  ~  l 

T  BAR  =  0«  0 

SUM-O. 0 

DO  1  I s  1 «  M 

SUM-SUM  +  XNl I  ) 

1  TBAR*XN( I  )*T  II  I ) ♦TBAR 
TBAR*T  BAR/ SUM 
RETURN 
END 

C  ZAP 

SUBROUTINE  ZAP<A,MAXZ) 

C  ZAP  RESCALES  DENSITIES  IN  THE  Z  DIRECTION. 

DIMENSION  A(1),TA«200) 

DO  44  K  =  1 , 4 

*  DO  40  J  =  1 »  MAXZ 
40  T  A (  Jl )  =  0  3  0 

Nl  —  C  MAXZ-3  )/Z 
L=N/2 ♦ 1 

DO  42  J-2f  MAXZ  «  2 
M»J*U-1)*MAXZ*1 
L  s  L  ♦  1 

42  T  A l L )  =  A l M ) 

DO  44  J si « MAXZ 
I s J ♦ l X" 1 ) ♦MAXZ 
44  A,  l  I  )  =  T  A  (  J  ) 

RETURN 

END 

*  LIST 

*  DATA 


PI* 

2  m  OOO 

P2  * 

Y  K  A  B  = 

2m 80E- 

30 

YK  B  A= 

YKBC* 

I  •  2  0  E  - 

29 

YKCB  = 

ZO* 

•  BOOO 

RO* 

TEMP* 

334.  0 

DNA  = 

1L.OOE 

12 

DNB* 

XK  OA  = 

3.67000 

XK  08  = 

E  = 

1L.  39E 

01 

EON* 

VA  = 

2.37E 

04 

VB= 

DAR  = 

4.91E 

01 

OBR* 

0  AZ  = 

4.  91 E 

01 

OiBZ* 

R  A0  = 

9.  3  8E 

03 

RBA  = 

R  BC  = 

4.  02 E 

04 

RCB* 

XNl* 

5.  79 E 

16 

XN2  = 
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0000 

2.00E- 

14 

1.68E- 

-13 

2500 

XL* 

7.50000 

1.0  OE 

00 

DNC* 

l.OOE 

00 

05000 

XKOC* 1 

.8  5000 

2.40E- 

•16 

1.1 9E 

04 

VC* 

1.19'E 

04 

2.4BE 

01 

DCR* 

2  •  4  B  E 

01 

2.48E 

01 

DCZ* 

2.48E 

01 

1.16E 

03 

9. 726 

03 

0.0  OE 

00 
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A 

TIME 

1.4626E- 

10 

2.8483E-05 

409  7E- 09 

2a  9749E -05 

1.2068E- 

■07 

3a 1015E-05 

2 a  5800E- 

•06 

3.2281E-05 

2.9160E- 

•05 

3.3547E-05 

3.4199E-04 

3. 481 3E -05 

2.4735E- 

•03 

3.6079E-05 

1.7308E-02 

3.7345E-05 

9* 4264E- 

•02 

3, 861 IE-05 

3  j  8 108E-Q1 

3.9877E-05 

i • 3249E 

00 

4.  U43E-05 

4.2595E 

00 

4.2409E-05 

1.1447E 

01 

4. 3674E-05 

2.8971E 

01 

4a  4940E-05 

6o  4955E 

01 

4o  6206E-05 

1  a  4465E 

02 

4.7472E-05 

5a  2214E 

02 

5a  0004E-05 

U4318E 

03 

5.2536E-05 

3a 1407E 

03 

5#  5068E-05 

5#  4105E 

03 

5a  7600E-05 
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00 

O.OOOOE 

00 

6. 9551E-22 

4. 9308E-04 

O.OOOOIE 

00 

O.OOOOE 

00 

O.OOOOE 

00 

O.OOOOE 

00 

3. 1745E-23 

5. 1333E-04 

O.OOOOIE 

00 

O.OOOOE 

00 

O.OOOOE 

00 

O.OOOOE 

00 

1.6573E-24 

5. 2346E-04 

O.OOOOE 

00 

O.OOOOE 

00 

O.OOOOE 

00 

O.OOOOE 

00 

7. 4927E-26 

5.3359E-04 

O.OOOOE 

00 

O.OOOOE 

00 

O.OOOOE 

00 

O.OOOOE 

00 

3.2178E- 27 

5. 5384E-04 

O.OOOOE 

00 

O.OOOOE 

00 

O.OOOOE 

00 

O.OOOOE 

00 

1. 7373E-28 

5. 6397E-04 

O.OOOOE 

00 

O.OOOOE 

00 

O.OOOOE 

00 

O.OOOOE 

00 

8. 0226E” 30 

5.7410E-04 

O.OOOOE 

00 

O.OOOOE 

00 

O.OOOOE 

00 

O.OOOOE 

00 

3. 2679E- 3 1 

5.  943  5E-04 

O.OOOOE 

00 

0. OOOOE 

00 

O.OOOOE 

00 

O.OOOOE 

00 

1. 8171E”32 

6.0448E-04 

SUM  A  =  5*  80*796-02! 

SUMB  =  2.7725E-02 

SUMC*  9.1 773E-02! 

T  A=  1.2771E-04  TB*  1.7493IE-04  TC*  1.9557E-04 
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519164 


T HEORET  IC AIL  CURVES  -  THE  Y  AXIS  APPLIES  TO  ION  SPECIES  A 
0.34E  04+  A  B 


C 


0  a  2 7  E  04 


0.20E  04* 


0.I4E  044 


0.68E  03 


A  (5  A 


C  b 

B 

C  B 

B  B 

B 


B  C 

A  C 

B  C 
B  C 


B  » 


A  A  A  A 

D-°”?£-03  0.12F-03  0.I3E-03  0.  ISE-03  0.  I6E-03  0.iah-03  0.20E-03  0„2iE-03  0.23E-03  0.24E-03  0.26E-03  0„27E-03 


0^29E-0 


3 
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1 


519165 


A 

TIME 

B 

5.33376-07 

l • 43686-04 

1. 33436-04 

9a  43406-05 

1.4875E-04 

6.56406-03 

6.31276-03 

1 • 5 381 E -04 

1.06476-01 

1  •  66561:-  01 

1 . 5887 E -04 

9.94076-01 

3. 125 IE  00 

1.6394E-04 

2.27746  01 

2.68436  01 

1.6900E-04 

5.72  49E  01 

1.6446E  02 

1.7407E-04 

9.7683E  01 

5.38216  02 

1.7913 E -04 

1.23336  02 

1.12516  03 

1.8419E-04 

1. 0951 E  02 

l  •  44026  03 

1 • 8926E-04 

9.52736  01 

1  .  12491:  03 

1. 9432E-04 

8.49366  01 

5*84521:  02 

1. 99386-04 

7.67096  01 

2.23366  02 

2. 0445E-04 

7®  65.29E  01 

7.4382EE  01 

2. 095  IE -04 

7.79446  01 

3.8734E  01 

2.1457E-04 

8.38716  01 

2.9229EE  01 

2. 19646-04 

8.7565E  01 

2*  6628 E  01 

2. 24706-04 

9.15186  01 

2 •  5 3 5 2 EE  OIL 

2.29776-04 

1.0006E  02 

2.47136  01 

2.3483E-04 

1.0459E  02 

2 *437 If:  01 

2.39896-04 

1.09296  02 

2*  420 SEE  01 

2.44966-04 

1.19336  02 

2  •  4 1 3  7  EE  01 

2.5002E-04 

1.2459E  02 

2  •  4109EE  01 

2.5508E-04 

1. 29976  02 

2*41 16EE  01 

2.60156-04 

1.41036  02 

2* 41L06EE  01 

2.70286-04 

1.46706  02 

2  •  3  8  6 1  EE  01 

2. 8040E-04 

1.52  51E  02 

2. 3320E  01 

2.90536-04 

1.65176  02 

2* 2461  EE  01 

3.00666-04 

1.72006  02 

2.1884EE  01 

3.05726-04 

1.78826  02 

2 •  1242EE  01 

3.10796-04 

1.89406  02 

1.9619EE  01 

3.20916-04 

1.90086  02 

1.7530EE  01 

3.31046-04 

1.86256  02 

1.49036  OIL 

3.41176-04 

1.76816  02 

1* 173 IE  01 

3.51296-04 

1.31996  02 

9.90236  00 

3. 5636E-04 

8.28696  01 

8 *2421  EE  00 

3. 6142E-04 

4.05936  01 

5.1367EE  00 

3.71556-04 

1.35776  01 

3.0290E  00 

3.8168E-04 

1.4269E  00 

1.2835EE  00 

3.91806-04 

2®  69 72 E -01 

5  •  8  6  5  2  EE  -  0 1 

4.01936-04 

3.6944E-02 

1  •  5  2 1 2  EE  -  0 1 

4. 1206E-04 

1.8947E-03 

4.42586-02 

4.2219E-04 

2.15546-04 

6.9&90E-03 

4.3231E-04 

1.95276-05 

2 .80226-03 

4. 4244E -04 

7.66966-07 

3* 176  IE- 04 

4.52576-04 

6.91026-08 

5.5605E-05 

4.6270E-04 

5.15016-09 

5 . 2  5  8  7  E  -  06 

4.72826-04 

1.80136-10 

2.2256E-06 

4.82956-04 

1.44516-11 

1.5895E-07 

4.93086-04 

9.68076-13 

2  *  2469E-08 

5. 0321 E -04 

3.18116-14 

1* 6462 E- 09 

5.1333E-04 

2.41386-15 

7. 440 IE -10 

5.23466-04 

1.52576-16 

O.OOOOE  00 

O.OOOOE  00 

4.79726-18 

O.OOOOE  00 

O.OOOOE  00 

3.5454E-19 

O.OOOOE  00 

O.OOOOE  00 

2.16556-20 

O.OOOOE  00 

O.OOOOE  00 

4.43406-21 

O.OOOOE  00 

O.OOOOE  00 

O.OOOOE  00 

O.OOOOE  00 

O.OOOOE  00 

O.OOOOE  00 

O.OOOOE  00 

O.OOOOE  00 

O.OOOOE  00 

TIME 

C 

TIME 

1.53816-04 

3.6475E- 

•05 

1. 538  IE-04 

1.58876-04 

1.64426- 

03 

1. 58876-04 

1.63946-04 

2.58386- 

■02 

1.63946-04 

1. 69006-04 

1.63606 

00 

1.7407E-04 

1.791L3E-04 

7.65666 

00 

1.79136-04 

1.84196-04 

2.3578E 

01 

1.84196-04 

1.89266-04 

5.3304E 

01 

1.89266-04 

1.9432E-04 

1.38846 

02 

1.99386-04 

2.04456-04 

1.76616 

02 

2. 0445E-04 

2.09516-04 

2.06496 

02 

2.09516-04 

2.14576-04 

2.48666 

02 

2.19646-04 

2.2470E-04 

2.66036 

02 

2.24706-04 

2.29776-04 

2. 82876 

02 

2.2977E-04 

2.34836-04 

3.18726 

02 

2.3989E-04 

2.44966-04 

3.37116 

02 

2.4496E-04 

2.5002E-04 

3.55636 

02 

2. 5002E-04 

2.55086-04 

3.93176 

02 

2.601 5E-04 

2.6521E-04 

4.12346 

02 

2.65216-04 

2.70286-04 

4.3191E 

02 

2.70286-04 

2.75346-04 

4.52206 

02 

2.7534E-04 

2.8547E-04 

4.7337E 

02 

2.80406-04 

2.9053E-04 

4.9516E 

02 

2.85476-04 

2.9559E-04 

5. 1754E 

02 

2.90536-04 

3.0572E-04 

5.4056E 

02 

2.95596-04 

3. 1079E-04 

5.64126 

02 

3.00666-04 

3.15856-04 

5.8791E 

02 

3.05726-04 

3.2598E-04 

6.12056 

02 

3.10796-04 

3.31046-04 

6.36806 

02 

3.15856-04 

3.36106-04 

6.63266 

02 

3.2091E-04 

3.4623E-04 

6.9126E 

02 

3.2598E-04 

3.51296-04 

7.2050E 

02 

3.  310  4E-04 

3.56366-04 

7.5053E 

02 

3.3610E-04 

3.61426-04 

7.7751E 

02 

3.411 7E-04 

3.7155E-04 

7.96336 

02 

3.462  36-04 

3.81686-04 

8.0254E 

02 

3.51296-04 

3.91806-04 

7.37056 

02 

3.61426-04 

4.01936-04 

5.57956 

02 

3.71556-04 

4.22196-04 

3.5167E 

02 

3.8168E-04 

4.32316-04 

1.6214E 

02 

3.91806-04 

4.4244E-04 

2.01036 

01 

4. 1206E-04 

4.62706-04 

4.63986 

00 

4.2219E-04 

4.7282E-04 

7.28196- 

01 

4.32316-04 

4.8295E-04 

3.4161E- 

02 

4.52576-04 

5.03216-04 

4.43386- 

03 

4. 6270E-04 

5.13336-04 

4.3194E- 

04 

4.72826-04 

5. 2  346E-04 

1.45156- 

05 

4. 9308E-04 

5.4372E-04 

1.47996- 

06 

5.03216-04 

5.53846-04 

1. 16596- 

07 

5.13336-04 

5.63976-04 

3.37906- 

09 

5. 33596-04 

5.8423E-04 

3.09896- 

10 

5.437 2 E- 04 

5.94356-04 

2. 19566- 

11 

50  5 3 84 E- 04 

6.04486-04 

5.8150E- 

13 

5. 7410E-04 

6.24746-04 

5.10896- 

14 

5.8423E-04 

6.34866-04 

3.42056- 

15 

5.94356-04 

6. 4499E-04 

2.0067E- 

16 

6.0448E-04 

6.5512E-04 

8.43676- 

17 

6.14616-04 

0.00006  00 

7. 3686E- 

18 

6.2474E-04 

O.OOOOE  00 

4.7880E- 

19 

6.  3486E-04 

O.OOOOE  00 

2.72186- 

20 

6. 4499E-04 

25  JUL  72  PAGE 


2 


519166 


Oo  OOOOE 

00 

0»  OOOOF 

00 

0®  0  000  E 

CIO 

0. OOOOE 

00 

1.1008E-20 

6*  55 1 2E -04 

0. OOOOE 

00 

0.  OOOOE 

00 

0. OOOOE 

00 

□•OOOOE 

00 

9.7908E-22 

6*  6524E -04 

SUM  A  3 

3.0290E-02 

SUMfl.* 

2*  62221b -02 

S  UMC  = 

9.7792E-02 

T  A=  1 *  9666E- 04 

TB  3 

2 o  94786- 04 

T  C  = 

3.0584E-04 

25  JUl  72  PAGE 


THEORETICAL.  CURVES  «  THE  Y  AXIS  APPLIES  TO  ION  SlPECIES  A 


0.14E  04* 


Oo  12E  04* 


0.86E  03+ 


0.58E  03 ♦ 


0«»  2  9E  03  + 


519167 


OiO0E  0()i _ * _ +  _^_ . + _ 4  4  4  '»  *  *  A  A  A  A  c 

0.1liE-l)3  0.18E-03  0.20E-CI3  0.22E-03  0.25E-03  0.27E-03  0..29IE-03  0.31E-03  0.34E-02I  0.36E-03  0. 38E-03  ~0o4ClE-03  0.<,2E-0 
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519168 


A 

TIME 

B 

TIME 

C 

TIME 

4.5954E- 

22 

1.7407E-04 

7*73276 

-23 

1*74076-04 

2*21916-19 

1*79136-04 

4.1456E- 

17 

1.7913E-04 

8*44236 

-16 

1*84196-04 

3*83346-16 

1*84196-04 

I.9250E- 

14 

1.8419E-04 

4*07786 

-13 

1*  8 9 26 E-04 

1*33506-13 

1*89266-04 

1*7574E- 

11 

1. 89266-04 

1 *  1 660 E 

-10 

1® 9432E-04 

1*56906-11 

1*  9432E  — 04 

1.5097E- 

09 

1 . 9432 E -04 

7 *  0647E 

-09 

1*99386-04 

8*71716-00 

2*  0445E-04 

2.3672E- 

07 

1.9938E-04 

1*21686 

-05 

2*09516-04 

3*09176-06 

2*09516-04 

8.1500E- 

06 

2.0445E-04 

3*  51 04E 

-04 

2* 1457E-04 

6* 1846E-05 

2*14576-04 

4.1711E- 

04 

2.0951E-04 

4*60646 

-03 

2*19646-04 

1*  3159E— 02 

2*24706-04 

7.3709E- 

03 

2.1457E-04 

3*97526 

-01 

2*2977E-04 

1*07506-01 

2*29776-04 

1.4920E- 

01 

2.1964E-04 

2* 1 707E 

00 

2*  348 3 E -04 

5*86906-01 

2*34836-04 

1.3845E 

00 

2.2470E-04 

7*61756 

00 

2*  3989E-04 

8*32436  00 

2*4496E-04 

1.1479E 

01 

2.2977E-04 

3*  8000E 

01 

2»5002E-04 

1*  98I56E  01 

2*50026-04 

5.3827E 

01 

2.3483E-04 

5*  3 1 64 E 

01 

2*55086-04 

3*72396  01 

2*55086-04 

1.8559E 

02 

2.3989E-04 

5*96  796 

01 

2*  60 15  E-04 

7*85346  01 

2*65216-04 

4.1131E 

02 

2.4496E-04 

5*  7034E 

01 

2*65216-04 

9*60006  01 

2*7028E-04 

6.0617E 

02 

2.5002E-04 

5* 1040E 

01 

2*  70 28 E-04 

1*10336  02 

2*75346-04 

6*  0802 E 

02 

2m 55086-04 

4.5669E 

01 

2*75346-04 

1*  3230E  02 

2*85476-04 

4  •  i  6  0  0  E 

02 

2.6015E-04 

4*  2 095E 

01 

2*80406-04 

1*41516  02 

2*90536-04 

8.5318E 

OIL 

2. 7028E-04 

4*  03 06 E 

01 

2® 8 5 47 E-04 

1*50176  02 

2*  9559E-04 

2.3320E 

OIL 

2.8040E-04 

3*  98196 

01 

2® 9053E-04 

1*67086  02 

3*05726-04 

1  •  8  1 8  0  E 

OIL 

2. 9053E-04 

4*  01 54E 

01 

2*95596-04 

1*75536  02 

3* 1079E-04 

1.7526E 

OIL 

3.0066E-04 

4. 1052E 

01 

3* 0066 E-04 

1*84006  02 

3* 1585E-04 

1  •  7  5  5  0  E 

OIL 

3.0572E-04 

4.2352E 

01 

3*05726-04 

2*01006  02 

3*25986-04 

1  •  7  6  8  9  E 

OIL 

3. 1079E-04 

4*  3939E 

01 

3*10796-04 

2*09496  02 

3*  3104E-04 

1.8148E 

OIL 

3. 2 09 IE -04 

4*56846 

01 

3*15856-04 

2*18096  02 

3«3610E -04 

1.873IE 

OIL 

3.31 04E -04 

4*75336 

01 

3*20916-04 

2*36036  02 

3*46236-04 

1.9313E 

OIL 

3.4117E-04 

4*  9445E 

01 

3*25986-04 

2*47696  02 

3*  51 2  9E -04 

1.9834E 

OIL 

3.5129E-04 

5® 1 409E 

Oi 

3*  3 104E-04 

2*59586  02 

3*  5636E-04 

2.0069E 

01 

3.5636E-04 

5  «  3  4  3  7  E 

01 

3*  3610  E-04 

3*01966  02 

3*  715 5E -04 

2.0267E 

01 

3. 6 142E-04 

5*  5  5  64 E 

01 

3*41176-04 

3* 653 IE  02 

3*91806-04 

2.0549E 

01 

3.7155E-04 

5*  7842E 

01 

3*  4623E-04 

3* 990 IE  02 

4*01936-04 

2 • 069  3 E 

OIL 

3.8168E-04 

6*  0451 E 

01 

3*51296-04 

4*  3606E  02 

4*12066-04 

2.0647E 

OIL 

3.9180E-04 

6*33166 

01 

3*56366-04 

4*74036  02 

4*22196-04 

2.0425 E 

OIL 

4.0193E-04 

6*  66  04 E 

01 

3*  61 42 E-04 

5*56966  02 

4»4244E-04 

1.9943E 

01 

4. 1206E-04 

7  #  3  3  8 1 E 

01 

3* 7 155 E-04 

6*  OO06E  02 

4*52576-04 

1.9264E 

01 

4.2219E-04 

8.0638E 

01 

3*8 168  E-04 

6*  4308E  02 

4*62706-04 

1.8412E 

OIL 

4.3231E-04 

9*  6447 E 

01 

4*01936-04 

6*  226 5 E  02 

4*  8295E  — 04 

1  a  7  3  7  5  E 

OIL 

4.4244E-04 

1*  0501 E 

02 

4* 1206 E-04 

5* 2276E  02 

4*  9308E  -04 

1.5B59E 

01 

4.5257E-04 

1*13946 

02 

4*22196-04 

3*  7 13 BE  02 

5*03216-04 

1.391  IE 

OIL 

4. 627 OE -04 

1*  33276 

02 

4*42446-04 

1*  0510E  02 

5*  2346E  — 04 

1.1487E 

OIL 

4.7282E-04 

U4357E 

02: 

4»5257E-04 

4* 1 173 E  01 

5*  3359E-04 

0.8O74E 

00 

4.8295E-04 

1  *  5  2  7  5  E 

02 

4* 62 70 E-04 

1*191 IE  01 

5* 43 7 2E— 04 

6.0783E 

00 

4. 93086-04 

1.4763E 

02 

4® 8295 E-04 

1*1424E  00 

5*  6397E  — 04 

3.7322E 

00 

5.0321 E -04 

1.2248E 

02 

4*  9308  E-04 

2*  4060E-0 1 

5*74106-04 

2.0093E 

00 

5.1 333E-04 

8  •  7  6  3  6  E 

01 

5*03216-04 

3*  92876—02 

5*  8423E-04 

9.0235 E- 

OIL 

5. 23466-04 

2.52 14E 

01 

5*23466-04 

2*24386-03 

6*04486-04 

3.8816 E- 

■OIL 

5.335 9E ~04 

1.0659E 

01 

5*  33 59 E-04 

3*  36 15 E-04 

6*14616-04 

1.0918E-01 

5.4372E-04 

3.42806 

oo 

5*43726-04 

4*02236-05 

6*  2474E-04 

4 . 1  7  5  3  E  - 

02 

5. 5384 E -04 

3*18646 

-01 

5*63976-04 

1*74496-06 

6«  4499E-04 

1.0271  E:- 

02 

5. 63976-04 

7.8127E 

-02 

5*74106-04 

2* 1526E-Q7 

6*55126-04 

4.212^E- 

04 

5. 8423E-04 

1.4379E 

-02 

5*  8423  E-04 

2*13846-08 

6*  6524E-04 

1.4727E- 

-04 

5 . 943  5E -04 

7.3326E 

-04 

6*04486-04 

3*76456-09 

6*  7537E-04 

2.3593E- 

05 

6.0448E-04 

1 .3126E 

-04 

6* 1461 E-04 

7*82486-10 

6*85506-04 

O.OOOOE 

00 

O.OOOOE  00 

1.7567E 

-05 

6*  2 4 74 E-04 

8*61626-11 

6*95636-04 

0  «  OO'OOE 

00 

O.OOOOE  00 

1.9585E 

-06 

6*  3486  E-04 

7*60936-12 

7*  0575E  -04 

O.OOOOE 

00 

O.OOOOE  00 

9.5632E 

-08 

6*  55 12 E-04 

1*06146-12 

7* 15 80 E-04 

0  a  0  0  0  0  E 

00 

O.OOOOE  00 

1.0575E 

-08 

6*  65 24 E-04 

2*  46896-  I  3 

7  »  2  6  0 1 E  -  04 

O.OOOOE 

00 

0®  OQOOE  00 

9. 8961 E 

-10 

6*  7537  E-04 

2*  5237E- 14 

7*  3614E-04 

O.OOOOE: 

00 

O.OOOOE  00 

3.05266 

-10 

6*85506-04 

0*  OOOOE  00 

O.OOOOE  00 

O.OOOOE  00 
O.OOOOE  00 
O.OOOOE  00 
O.OOOOE  00 

SUMA-  1.7250E-02 
SUMB=  2. 3678E-02 
SUHC=  9.2872E-02 

T  A  =  2.7064E-04  TB«  21.81L45E-04 


O.OOOOE  00  4.2053E 

O.OOOOE  00  4.1291E 

O.OOOOE  00  3.  4415  5  E 

O.OOOOE  00  9.98I33E 


Cn 

C* 


25  JUL  72  PAGE 


2 


519169 


11 

6.9563E-04 

O.OOOOE 

00 

O.OOOOE 

00 

12 

7.0575E-04 

O.OOOOE 

00 

O.OOOOE 

00 

13 

7. 1588E-04 

O.OOOOE 

00 

0.  OOOOE 

00 

14 

7.2601 E-04 

O.OOOOE 

00 

0. OOOOE 

00 

TO  4.0497E-04 


25  JUL  72  PAGE 


1 


519170 


THEORETICAL  CURVES  -  THE  Y  AXIS  APPLIES  TO  ION  SPECIES  A 
0 <i  6 i E  03+  A  A 


0.49E  03* 


B 

C 


c; 

B 


D 


0.36E  03 ♦ 


£  0  a  2  4  E  03^ 


0.12E  03  ♦ 


A 


A 


D 


B 

C 


C 


B 


B 

B  DC 


13  B  B 
C 


80 
D  C 


B 

BC 

C 


B 

C 


C 


c 


D 


D 


I  A 

I  C 

i  a.  c; 

I  ABC 

J,  03E  OO* - ♦ - ♦ 

Oj21r - J  3  0#  2  4E-03  0  a  2  7  F 


A  A  A  AA  A 

A  A  A  A  A 

- ♦ - ♦ - ► - 

03  0,  jOE-C3  0.J3E-03  0ojblE-03 


A  A  A  A  A  A  A 

A  A 

- 4 - - - 4 - - 4 - 

0.39E-03  0,42E-03  0.44F-03 


D 

AAA  D 

- 4 - + - 4 - 

0.47E-03  0 , 50E- 03  0.33E-03 


0..56IE-0 


3 


25  JUL  72  PAGE 


1 


519171 


A 

TIME 

0 

2  •  09fi»3E 

34 

2.0445E-04 

1.3376E- 

”06 

W8809E 

--30 

2a  095  IE -04 

2.4874E 

-04 

8.3094E 

-27 

2.1457E-04 

3.0048E 

-03 

8.5566E 

-24 

2* 1964E-04 

2.2078E-02 

5*7  42!  OE 

-21 

2.2470E-04 

7.728QE-01 

1.9145E 

-18 

2#  2977E-04 

2.767  9E 

00 

4.0292E 

“16 

2* 3403E -04 

7.235  4E 

00 

5.9632E 

-14 

2.3989E-04 

2  a  2  4  9 1 E 

01 

5.  3172E 

”12 

2a4496E-04 

2a 763 OE 

01 

4#  040 IE 

”10 

?a  5  002E--04 

2a  8720E 

01 

!• 8025E 

”08 

2#  550  8E-04 

2a  42  37E 

01 

7*  5789E 

”07 

2  a  6  0 1 5  E  -  04 

2.2137E 

01 

la  8549E 

-05 

2. 652  IE ”04 

2.0857E 

01 

4.4O04E 

-04 

2#  7 02 8E” 04 

2.0304E 

01 

8.1141E 

-03 

2.7534E-04 

2a  0693E 

01 

8  a  0 1 5  3  E 

“02 

2.8040E-04 

2.3997E 

01 

6.2826E 

”01 

2a 8 547E -04 

2  •  6  2  8  3  E 

01 

4.3J78E 

00 

2#  9053E-04 

2a  88 19E 

01 

1.8790E 

01 

2a  9559E -04 

3a 1 572E 

01 

6.4803E 

01 

3.0066E-04 

3a  7925E 

01 

L.4943E 

02 

3a  0572E-04 

4.1336E 

01 

2#  492  IE 

02 

3.1079E-04 

4.4855E 

01 

2  a  9  6 1 1 E 

02 

3a 1505E-O4 

5.2839E 

01 

2  a  4  7 4  9  E 

02 

3* 2 09 IE ”04 

5a  7686 E 

01 

1.5574E 

02 

3a  2598E”04 

6.3167E 

01 

7 a  7942E 

01 

3.3104E-04 

7.5872E 

01 

3.4855E 

01 

3.3610E-04 

8.2708E 

01 

la  9166E 

01 

3#  41 17E-04 

8.9845E 

01 

la  3  6  9  9  E 

01 

3.4623E-04 

la  0537E 

02 

1.2212E 

01 

3. 5129E-04 

1  a  1 3  5  8  E 

02 

l a  158  2 E 

01 

3a  5636E-04 

1.2186E 

02 

1.1573E 

01 

3.  6142E-  04 

1.3277E 

02 

l a  173  IE 

01 

3.71556-04 

1.2789E 

02 

1.2167E 

01 

3a  8 168E-04 

1.1353E 

02 

1.2729E 

01 

3*9180E“04 

5.7701E 

01 

1.3334E 

01 

4a  0193E-04 

3.3312E 

01 

1.3940E 

01 

4* 1206E-04 

1.5665E 

01 

U4532E 

01 

4a  221 9E- 04 

2.6708E 

00 

1.5L02E 

01 

4a  323 IE” 04 

8a  9040E- 

”01 

l  a  5  6  4  3  E 

01 

4.4244E-04 

2.3505E-01 

1.6147E 

01 

4a  5257E-04 

2.02 03E-02 

1  a  6  5 1 1 E 

01 

4. 627 OE” 04 

4.6605E* 

“03 

1  a  7  0 1  0  E 

01 

4.7282E-04 

8.5023E-04 

1.7377E 

01 

4a8295E-04 

4. 8567E-05 

W7697E 

01 

4.9308E-04 

8.09O4E-O6 

1  a  7  9  7  7  E 

01 

5, 0  32  IE -04 

1.2750E- 

-06 

U8098E 

01 

5.  1333E-04 

O.OOOOE 

00 

1  a  8  0  3  6  E 

01 

5a  2 346E -04 

OaOOOOE 

00 

1.7759E 

01 

5.3359E-04 

OaOOOOE 

00 

1.7292E 

01 

5,  437 2E  -  04 

O.OOOOE 

00 

W6565E 

01 

5.5384E-04 

OaOOOOE 

00 

l *  5596E 

01 

5a6397E-04 

O.OOOOE 

00 

1.4296E 

01 

5.  741 OE- 04 

O.OOOOE 

00 

la  2699E 

01 

5.  8423E-04 

0.00006 

00 

1.0748E 

01 

5.9435E-04 

0. QOQOE 

00 

8  a  5  3  4  9  E 

00 

6.0448E-04 

O.OOOOE 

00 

6.2848E 

00 

6a 1 46 1 E - 04 

O.OOOOE 

00 

40 1640E 

00 

6.2474E-04 

O.OOOOE 

00 

2a  4042E 

00 

6.3486E-04 

O.OOOOE 

00 

TIME 

c 

TIMIE 

2.6521E-04 

6.7209E- 

■05 

2.7534E-04 

2. 7534E-04 

6. 0 101 E—04 

2#  8040E— 04 

2. 8040E--04 

4.  283 IE- 

■02 

2.9053E-04 

2.8547E--04 

2.3659E-  01 

2.9559E-04 

2. 9559E--04 

9.2925E-01 

3.0066E-04 

3a  0066E-04 

7.3203E 

00 

3. 1 079E-04 

3a  0572E-04 

1.4550E 

01 

3.  1585E-04 

3a 1585  E—04 

2. 3912E 

01 

3.2091E-04 

3*  209  IE-04 

4.3371E 

01 

3.3104E-04 

3, 2598  E-04 

5. 1305E 

01 

3.3610E-04 

3.3610E-04 

5. 7905E 

01 

3.4U7E-04 

3.4117E-04 

6.8456E 

01 

3.  5129E-  04 

3. 4623E-04 

7.3063E 

01 

3.  563 6 E-  04 

3a5636E-04 

7.7525E 

01 

3.6142E-04 

3.61426-04 

9.5693E 

01 

3.8168E-04 

3.8168E-04 

1.0563E 

02 

3a 91806- 04 

3.9180E-04 

1.1627E 

02 

4#  0193E-04 

4*  0 193 E-04 

1.4072E 

02 

4.2219E-04 

4a 1206E-04 

la  5430E 

02 

4. 32 31 E-04 

4* 323  IE- 04 

1  • 6  8  51 E 

02 

4.4244E-04 

4.4244E-04 

1.9894E 

02 

4. 6270E— 04 

4, 5257E—04 

2.1658E 

02 

4, 7282E- 04 

4*  7282  E—04 

2. 3666E 

02 

4.8295E-04 

4a  8295E-04 

2.8617E 

02 

5.  032  IE- 04 

4. 9308E-04 

3a 1370E 

02 

5. 1333E-04 

5* 1333E-04 

3.4252E 

02 

5. 2346E- 04 

5*  2 346 E-04 

4.0536E 

02 

5.4372E-04 

5.3359E-04 

4.3917E 

02 

5.5384E-04 

5.5384E-04 

4.7415E 

02 

5. 6397E-04 

5.6397E-04 

5.4295E 

02 

5. 8423E- 04 

5.7410E-04 

5.5785E 

02 

5.9435E-04 

5. 9435 E-04 

5.4213E 

02 

6. 0448E- 04 

6*0448E»04 

4.8566E 

02 

6.1461E-04 

6.1461 E-04 

2. 4440 E 

02 

6. 3486E-04 

6a  3486E— 04 

1.3764E 

02 

6.4499E-04 

6. 449 9 E-04 

6. 1116E 

01 

6. 5512E—04 

6.5512E-04 

1. 044 2E 

01 

6.753  7E-04 

6a  7537  E—04 

3« 163  IE 

00 

6.8550E-04 

6.8550E-04 

7.6217E- 

■01 

6.9563E-04 

6.9563E-04 

6*9304E- 

■02 

7.  1588E-04 

7.1508E-O4 

1. 408  IE- 

■02 

7.2601E-04 

7a  260 l E-04 

2. 3397E-03 

7.3614E-04 

7.3614E-04 

la  4880E- 

04 

7.5639E-04 

7.5639E-04 

2a  3586E-05 

7.6652E-04 

7. 6652E-04 

3.0915E-06 

7.7665E-04 

7. 7665E-04 

1.5635E-07 

7.9690E-04 

O.OOOOE  00 

O.OOOOE 

00 

O.OOOOE  00 

O.OOOOE  00 

O.OOOOE 

00 

O.OOOOE  00 

O.OOOOE  00 

O.OOOOE 

00 

O.OOOOE  00 

O.OOOOE  00 

O.OOOOE 

00 

O.OOOOE  00 

O.OOOOE  00 

OaOOOOE 

00 

O.OOOOE  00 

O.OOOOE  00 

O.OOOOE 

00 

O.OOOOE  00 

O.OOOOE  00 

O.OOOOE 

00 

O.OOOOE  00 

OaOOOOE  00 

O.OOOOE 

00 

O.OOOOE  00 

OaOOOOE  00 

O.OOOOE 

00 

O.OOOOE  00 

O.OOOOE  00 

O.OOOOE 

00 

O.OOOOE  00 

O.OOOOE  00 

O.OOOOE 

00 

O.OOOOE  00 

O.OOOOE  00 

O.OOOOE 

00 

O.OOOOE  00 

OaOOOOE  00 

O.OOOOE 

00 

O.OOOOE  00 

25  JUL  72  PAGE  2 


I  •  3  050E  00  6#  44991:  -04 

0.0000E:  00 

O.OOOOE 

00 

0«>  OOOOE 

00 

0. OOOOE 

00 

6.0395E-01  6.  5 5 121: -04 

o.ooooe:  oo 

o.ooooe 

00 

O.OOOOE 

00 

0. OOOOE 

00 

2.4458E-01  <>.6!>24E-04 

O.OOOOE:  00 

o.ooooe 

00 

O.OOOOE 

00 

0.  OOOOE 

00 

9.2420E-02  6.7537IE-04 

o.ooooe:  oo 

o.ooooe 

00 

O.OOOOE 

00 

0. OOOOE 

00 

3.3573E-02  6.8550IE-04 

o.ooooe:  oo 

o.ooooe 

GO 

O.OOOOE 

00 

0. OOOOE 

00 

I3.9644E-G3  <>•  9563IE-04 

o.ooooe:  oo 

o.ooooe 

00 

O.OOOOE 

00 

0. OOOOE 

00 

2.5  379E-03  7*  0 1575 IE  ”04 

o.ooooe:  oo 

O.OOOOE 

00 

O.OOOOE 

00 

0. OOOOE 

00 

ib.6986E-04  7.  1  50  8  IE -04 

o.ooooe:  oo 

O.OOOOE 

00 

O.OOOOE 

00 

0. OOOOE 

00 

2.0068E-04  7.2601IE-04 

o.ooooe:  oo 

o.ooooe 

00 

O.OOOOE 

00 

0. OOOOE 

00 

3.  3  195E-05  7.3614IE-04 

o.ooooe:  oo 

o.ooooe 

00 

O.OOOOE 

00 

0. OOOOE 

00 

7.8159E-06  7  •  462  6IE  -04 

o.ooooe:  oo 

o.ooooe 

00 

O.OOOOE 

00 

0. OOOOE 

00 

SUMA-  1.0  3  01  IE -02 

SUMb«  2.1 448IE-02 

SU*C==  8.  6 2 25 IE -02 

T  A  =  3.  559  3E-04  TB  = 

5.1706E-04  TC= 

5.3179E-04 

519172 


25  JUL  72  PAGE 


519173 


THEORETICAL  CURVES  -  THE  Y  AXIS  APPLIES  TO  ION  SPECIES  A 
0.30E  034  A 


0*24E  03  + 


A  A 


D 

C  C 
B 

B 

C 

B  B 

C 

B 

C 


Ot 1 B  E  03* 


D 


B 

C 


A 


A 


C 


B 


cn 


0  » 1 2  E  03  + 


0  »  5  9  E  02 


B 

B  C 
D  C 


BB  B  Hi 
C  C 

A  C 

c  c 


A  B  A 

C  A  AAA  A  A  A  A  A  A 

I  A  B  C  C 

Oo  21E-33  + - + - ♦ - ♦ - 4 - 

0  3  2  7 1:  -  0  3  0.31E-J3  0,34E-03  0,  38E-03  0,41E-03 


AAAAAAAAAAAAAA 


- — - - + - —  - - *- - + - 

0»  45  E-03  0.48E--03  0,*52E-03  0.55E-03 


D 


B 

C 


B 

C 


AAA 

AAA  D 

- iy - ♦ - ♦ - 

0* 59E-03  0*  62E- 03  0.66E-03 


0, 69E-0 


3 


DISTRIBUTION  LIST 


No,  of  No*  of 

Copies  Organization  Copies  Organization 


2  Commander 

Defense  Documentation  Center 
ATTN:  TIPCR 
Cameron  Station 
Alexandria,  Virginia  22314 

4  Director 

Advanced  Research  Projects 
Agency 

ATTN:  Dr.  S.  Ruby 
Mr.  J.  Wade 
Dr.  D.  Mann 
Rdr  §  Opt  Tech 
Department  of  Defense 
Washington,  D.C.  20305 

1  Director 

Weapons  Systems  Evaluation 
Group 

ATTN:  Dr.  G.  Rosen 
Washington,  D.C.  20305 

1  Director 

Institute  for  Defense 
Analyses 

ATTN:  Dr.  E.  Bauer 
400  Army-Navy  Drive 
Arlington,  Virginia  22202 

6  Director 

Defense  Nuclear  Agency 
ATTN:  STRA  (RAAE) 
Washington,  D.C.  20305 

4  Director 

Defense  Nuclear  Agency 
ATTN:  STAP  (APTL)  -  2  cys 
STRA  (RAEV) 

DDST 

Washington,  D.C.  20305 


4  DNA  Information  and  Analysis 
Center 

TEMPO,  General  Electric  Company 
ATTN:  Mr.  D.  Archer 
Mr.  W.  Knapp 
Dr,  W.  Chan 
Mr.  A.  Feryok 
816  State  Street 
Santa  Barbara,  California 
93102 

3  Director 

Defense  Communications  Agency 
ATTN:  Code  H820,  Mr.  C.  Herr 
Code  540,  Mr.  J.  Krcek 
Code  340,  Mr.  W.  Dix 
Washington,  D.C.  20305 

1  Director- 

Defense  Intelligence  Agency 
Washington,  D.C,  20301 

1  Commander 

U.S.  Army  Materiel  Command 
ATTN :  AMCDL 
Washington,  D.C.  20315 

1  Commander 

U.S.  Army  Materiel  Command 
ATTN:  AMCRD,  Dr.  J . V . R . Kaufman 

Washington,  D.C.  20315 

1  Commander 

U.S.  Army  Materiel  Command 
ATTN:  AMCRD -T ,  Dr,  H ,M,  El-Bisi 
Washington,  D.C.  20315 

1  Commander 

U.S.  Army  Materiel  Command 
ATTN:  AMCRD -TV,  Mr.  H.  Cohen 
Washington,  D.C.  20315 


59 


DISTRIBUTION  LI  S  i 


No.  of  No.  of 

Copies  Organization  Copies  Organization 


1  Commander 

U.S.  Army  Materiel  Command 
ATTN:  AMCRD-TE 
Washington,  D.C.  203 1 S 

1  Commander 

U.S.  Army  Materiel  Command 
ATTN:  AMCRD-TP 
Washington,  D.C.  20315 

1  Commander 

U.S.  Army  Materiel  Command 
ATTN:  AMCRD-WN, 

Mr.  J.  Corrigan 
Washington,  D.C.  20315 

1  Commander 

U.S.  Army  Aviation  Systems 
Command 

ATTN:  AMSAV-E 

12th  and  Spruce  Streets 

St.  Louis,  Missouri  63106 

1  Director 

IJ.S.  Army  Air  Mobility  Research 
and  Deve 1 opment  Laborat  ory 
Ames  Research  Center 
Moffett  Field,  California 
94035 

3  Commander 

U.S.  Army  Electronics  Command 
ATTN:  'Lech  Docu  Cen 
AMSEL-RD 
Dr.  ii.  Bomke 

Fort  Monmouth,  New  .Jersey 
0  7703 

2  Director 

Atmospheric  Sciences  Laboratory 
ATTN:  Mr.  11.  Ballard 
Mr.  M.  Diamond 
White  Sands  Missile  Range 
New  Mexico  88002 


2  Commander 

U.S.  A  rmy  Missile  Comman d 
ATTN :  AMSM1-R 

Vi.  MI -RBI. 

k i  .L-  L.oac  A  r s e n a  1  ,  A 1  ab a ma 

“i i  8 1  1 6 

1  Commando r 

U.S .  Arno  Tank -Automot i ve 
Command 

ATTN :  AMSTA-RHFL 
Warren,  Michigan  48090 

2  Commander 

U.S.  Army  Mobility  Equipment 
Research  q  Development  Center 
ATTN:  Tech  Docu  Cen,  Bldg.  315 

AMSME-RZT 

fort  Bel  voir,  Virginia  22060 

1  Commander 

U.S.  A  r  '  !i  <n  i  t  ions  Command 

ATI  N  :  AMSMi 1  -  RE 

D" ver,  !T  v  Jersey  07801 

1  PIASTRE 

D.;,.  Arm}  Dicatinny  Arsenal 
ATTN :  SMUPA-FR-M-D 

i >  o  v e  r ,  N e w  <  J  e  r s  e y  0 7801 

1  (Commander 

U.S.  Anne  Picatinny  Arsenal 

AT  ITT  :  SMUPA-TW 

Dover,  .N c  v  Terse y  u 7 8 1 )  1 

2  Com;  and  or 

U.S.  .Trail}  Weapons  Command 
ATTN  :  AM  AWE -RE 
AMS WE -RDF 

Rock  Island,  Illinois  61202 


60 


DISTRIBUTION  LIST 


No.  of 

Copies  Organization 

2  Commander 

U.S.  Army  SAFEGUARD  System 
Command 
ATTN:  SSC-D 

SSC-HS,  Dr.  Ira  Merritt 
P.  0.  Box  1500 
Huntsville,  Alabama  35807 

2  Commander 

U.S.  Army  SAFEGUARD  System 
Office 

ATTN:  Mr.  J.  J.  Shea 
Mr.  J.  Davidson 
Commonwealth  Building 
1320  Wilson  Boulevard 
Arlington,  Virginia  22209 

1  Commander 

U.S.  Army  SAFEGUARD  System 
Command  Field  Office 
Bell  Telephone  Laboratories 
ATTN:  SSC-DEF-B,  Mr.  B. Turner 
Whippany  Road 

Whippany,  New  Jersey  07981 
1  Commander 

U.S.  Army  SAFEGUARD  System 
Evaluation  Agency 
ATTN:  EAB,  R.  DeKinder,  Jr. 
White  Sands  Missile  Range 
New  Mexico  88002 

1  Director 

U.S.  Army  Advanced  Materiel 
Concepts  Agency 
2461  Eisenhower  Avenue 
Alexandria,  Virginia  22314 

3  Commander 

U.S.  Army  Harry  Diamond 
Laboratories 
ATTN:  AMXDO-TD/002 
AMXDO-TI/012 
AMXDO-NP,  F.  Wimenitz 
Washington,  D.C.  20438 


No.  of 

Copies  Organization 

1  Commander 

U.S.  Army  Materials  and 
Mechanics  Research  Center 
ATTN :  AMXMR-ATL 
Watertown,  Massachusetts 
02172 

1  Commander 

U.S.  Army  Natick  Laboratories 
ATTN:  AMXRE,  Dr.  D.  Sieling 
Natick,  Massachusetts  01762 

1  Commander 

U.S.  Army  Satellite  Communi¬ 
cations  Agency 
ATTN:  Sys  Eng  Off 
Fort  Monmouth,  New  Jersey 
07703 

1  Chief 

U.S.  Army  Communications 
Systems  Agency 
Fort  Monmouth,  New  Jersey 
07703 

1  Commander 

U.S.  Army  Combat  Developments 
Command 

Nuclear  Agency 
ATTN :  Lib 

Fort  Bliss,  Texas  79916 
1  Commander 

U.S.  Army  Air  Defense  School 
ATTN:  Cmd  §  Staff  Dept 
Fort  Bliss,  Texas  79916 

1  HQDA  (DARD-DD) 

Washington,  D.C.  20310 

1  HQDA  (DARD-MSC-S) 

Washington,  D.C.  20310 


61 


DISTRIBUTION  LIST 


No.  of 

Copies  Organization 

1  HQDA  (DARD-ARE/CPT  K.  Miller, 
Dr.  F.  dePercin , 

Dr,  L,  Alpert) 
Washington,  D.C.  20310 

1  HQDA  (DACE) 

Washington,  D.C.  20310 

3  Commander 

U.S,  Army  Research  Office 
(Durham) 

ATTN:  Dr.  H.  Robl 
Dr.  A.  Dodd 
Dr.  R.  Mace 
Box  CM,  Duke  Station 
Durham,  North  Carolina  27706 

2  Director 

U.S.  Army  Advanced  Ballistic 
Missile  Defense  Agency 
ATTN:  Dr.  R.  Ruffine 

Dr.  J.  Gilstein 
1320  Wilson  Boulevard 
Arlington,  Virginia  22209 

2  Director 

U.S.  Army  Advanced  Ballistic 
Missile  Defense  Agency 
ATTN:  CRDABH-S,  Mr.  M.  Capp 
Mr.  H.  Solomons on 
P.  0.  Box  1500 
Huntsville,  Alabama  35809 

1  Chief  of  Naval  Operations 
ATTN:  OP -75 
Department  of  the  Navy 
Washington,  D.C.  20350 

3  Commander 

U.S.  Naval  Air  Systems 
Command 

ATTN:  AIR-604 
Washington,  D.C.  20360 


No.  of 

Copies  Organization 

3  Commander 

U.S.  Naval  Ordnance  Systems 
Command 

ATTN:  ORD-0632 
ORD-035 
ORD-5524 

Washington,  D.C.  20360 

3  Chief  of  Naval  Research 

ATTN:  Code  418,  Dr.  T.  Quinn 
Mr.  R.  Joiner 
Dr.  J.  Dardis 
Department  of  the  Navy 
Washington,  D.C.  20360 

2  Commander 

U.S.  Naval  Electronics 
Laboratory 
ATTN:  Mr.  W.  Moles 
San  Diego,  California  92152 

1  Commander 

U.S.  Naval  Ordnance  Laboratory 
ATTN:  Code  753,  Tech  Lib 
Silver  Spring,  Maryland  20910 

3  Director 

U.S.  Naval  Research  Laboratory 
ATTN:  Code  2020,  Tech  Lib 

Code  7700,  Dr.  K.  Hain 
Dr.  C.  Johnson 
Washington,  D.C.  20390 

2  Superintendent 

U.S.  Naval  Postgraduate  School 
Physics  Department 
ATTN:  Mr.  E.  Milne 
Dr.  0.  Heinz 

Monterey,  California  93940 

3  HQ  USAF  (AFNIN ;  AFRD ;  AFRDQ) 
Washington,  D.C.  20330 


62 


DISTRIBUTION  LIST 


No .  of 

Copies  Organization 

1  HQ  USAF  (AFTAC,  COL  C. 

Anderson) 

Washington,  D.C.  20330 

1  AFOSR  (SREP) 

1400  Wilson  Boulevard 
Arlington,  Virginia  22209 

2  AFSC  (SCTSW/LTC  R.  Linkous ; 

SCS) 

Andrews  AFB 

Washington,  D.C.  20331 

2  AFWL  (SYT,  Maj  W.  Whitaker 
1  LT  S.  Rockwood) 
Kirtland  AFB 
New  Mexico  87117 

5  AFCRL  (LKD,  Dr.  R.  Narcisi 
LKS,  Dr.  R.  Huffman 
DPR,  Dr.  A.  Stair 
LIJ,  Jr.  J.  Ulwick 
LI,  Dr.  W.  Pfister) 

L.  G.  Hanscom  Field 
Bedford,  Massachusetts  01730 

4  AFCRL  (LKB,  Dr.  K.  Champion 
Dr.  T.  Keneshea 
Dr.  J.  Paulson 
Dr.  E.  Murad 
L.  G.  Hanscom  Field 
Bedford,  Massachusetts  01730 

1  ARL  (ARC,  Dr.  T.  Tieman) 
Wright -Patterson  AFB 
Ohio  45433 

2  Director 
Institute  for 

Telecommunication  Sciences 
ATTN:  R44 ,  Mr.  G.  Falcon 
Mr.  G.  Adams 

Boulder,  Colorado  80302 


No .  of 

Copies  Organization 

1  Director 

National  Oceanic  and 

Atmospheric  Administration 
ATTN:  Code  179 

Mr.  T.  Ashenfelter 
U.S.  Department  of  Commerce 
Washington,  D.C.  20234 

5  Director 

National  Oceanic  and  Atmospheric 
Administration 
ATTN:  Dr.  E.  Ferguson 

Dr.  F.  Fehsenfeld 
Dr.  D.  Albritton 
Dr.  A.  Schmeltekopf 
Dr.  G.  Reid 

U.S.  Department  of  Commerce 
Boulder,  Colorado  80302 

5  Director 

National  Bureau  of  Standards 
ATTN :  Mr .  J .  Herron 
Mr.  D.  Garvin 
Mr.  B.  Steiner 
Dr.  H.  Weisman 
Reports  Section 
U.S.  Department  of  Commerce 
Washington,  D.C.  20234 

1  Director 

Argonne  National  Laboratory 
Library  Services  Department 
Report  Section,  Bldg  203 
RM-CE-125 

9700  South  Case  Avenue 
Argonne,  Illinois  60440 

1  Director 

Brookhaven  National  Laboratory 
ATTN:  Clas  Docu  Sec 
25  Brookhaven  Avenue 
Upton,  New  York  11973 


63 


DISTRIBUTION  LIST 


No .  of 

Copies  Organization 

1  Director 

Lawrence  Livermore  Laboratory 
ATTN:  Dr*  K.  Watson 
P.  0.  Box  912 

Berkeley,  California  94704 

3  Director 

Los  Alamos  Scientific 
Laboratory 

ATTN:  Library,  Dr*  H.Hoerlin 
Dr.  W.  Myers 
P.  0.  Box  1663 

Los  Alamos,  New  Mexico  87544 

4  Headquarters 

National  Aeronautics  and  Space 
Admini s t  rat i on 

ATTN:  Code  SG,  Dr.  M.  Dub in 
Code  SL,  Dr.  R. 

Horowitz 
Dr.  R.  Fellows 
Code  RRP ,  Dr.  A.  Gessow 
Washington,  D.C.  20546 

3  Director 

National  Aeronautics  and 
Space  Administration 
Ames  Research  Center 
ATTN:  Dr.  M.  Loewenstein 

Dr.  I.  Poppoff 

Dr.  R.  Whitten 

Moffett  Field,  California 
94035 

2  Director 

National  Aeronautics  and 
Space  Administration 
Goddard  Space  Flight  Center 
ATTN:  Code  691,  Dr.  M.  Mumma 

Code  625,  Dr.  Aiken 
Greenbelt,  Maryland  20771 


No .  of 

Copies  Organization 

1  Director 

National  Aeronautics  and 
Space  Administration 
Manned  Spacecraft  Center 
ATTN:  Dr.  R.  Hudson 

Houston,  Texas  77058 

3  Director 

National  Science  Foundation 
Atmospheric  Sciences  Section 
ATTN:  Dr.  E.  Droessler 

Mr.  L.  Craine 
Dr.  R.  Schoen 
1800  G  Street,  NW 
Washington,  D.C.  20550 

1  Aerochem  Research  Labs,  Inc, 

ATTN:  Dr.  H.  F.  Cal  cote 

P.  0.  Box  12 

Princeton,  New  Jersey  08540 

2  Beil  Telephone  Laboratories,  Inc. 
ATTN:  Tech  Rpt  Cen 

Dr.  L.  McKnight 
Whippany  Road 

Whippany,  New  Jersey  07981 

2  The  Dewey  Electronics  Corporation 
ATTN:  Dr.  M.  Hirsh 

Dr.  P.  Eisner 
11  Park  Place 

Paramus,  New  Jersey  07652 

1  Franklin  GNO  Corporation 
ATTN:  Dr.  M.  Cohen 
P.0.  Box  3250 

West  Palm  Beach,  Florida  33402 

2  GCA  Corporation 

ATTN :  Dr .  F .  Marmo 

Dr.  A.  Dalgarno 
Burlington  Road 
Bedford,  Massachusetts  01730 


64 


DISTRIBUTION  LIST 


No.  of 

Copies  Organization 

2  General  Electric  Company 
Valley  Forge  Space  Technology 

Center 

ATTN:  Dr.  M.  Bortner 

Dr.  T.  Baurer 
P.  0.  Box  8555 
Philadelphia,  Pennsylvania 
19101 

1  General  Research  Corporation 
ATTN:  Dr.  J.  Ise 
P.  0.  Box  3587 
Santa  Barbara,  California 
93105 

3  Gulf  General  Atomic,  Inc. 

ATTN :  Lib 

Dr.  V.  Van  Lint 
Dr.  B.  Turner 
P.  0.  Box  1111 

San  Diego,  California  92112 

1  Kitt  Peak  National  Observatory 
ATTN :  Library 
P.  0.  Box  4310 
Tucson,  Arizona  85717 

1  KMS  Industries 

ATTN:  Dr.  K.  Brueckner 
11689  Sorrento  Valley  Road 
San  Diego,  California  92121 

1  KMS  Industries 

ATTN:  Dr.  W.  Hooker 

1810  Burnet  Avenue 

Van  Nuys,  California  91405 

1  Lockheed  Missile  and  Space 
Company 

ATTN:  Dr.  L.  Fisher 
1111  Lockheed  Way 
Sunnyvale,  California  94086 


No.  of 

Copies  Organization 

4  Lockheed  Palo  Alto  Research 
Laboratory 
ATTN:  Library 

Dr.  J.  Kumer 
Mr.  R.  Sears 
Dr.  R.  Gunton 
3251  Hanover  Street 
Palo  Alto,  California  94304 

4  Mission  Research  Corporation 
ATTN:  Dr.  C.  Longmire 
Mr.  R.  Hendrick 
Mr.  P.  Fisher 
Mr.  M.  Schiebe 
812  Anacapa  Street,  Studio  5 
Santa  Barbara,  California 
93101 

1  Mitre  Corporation 
ATTN:  Tech  Lib 
P.  0.  Box  208 

Bedford,  Massachusetts  01730 

1  NRA,  Inc. 

ATTN:  Mr.  J.  Baer 
12  Nevada  Drive 

New  Hyde  Park,  New  York  11040 

3  The  Rand  Corporation 
ATTN:  Tech  Lib 

Dr.  C.  Crain 
Mr.  P.  Tamarkin 
1700  Main  Street 
Santa  Monica,  California 
90406 

2  R$D  Associates 

ATTN:  Dr.  F.  Gilmore 
Dr.  R.  LeLevier 
P.  0.  Box  3580 
Santa  Monica,  California 
90403 


65 


DISTRIBUTION  LIST 


No.  of  No.  of 

Copies  Organization  Copies  Organization 


4  Sandia  Corporation 

ATTN:  Org  3141,  Tech  Lib 
Org  1424 
Org  1433 

Org  100,  Mr.  Frank 
Hudson 

P.  0.  Box  5800 

Albuquerque,  New  Mexico  87115 

1  Science  Applications,  Inc. 

ATTN:  Mr.  R.  Lowen 

1250  Prospect  Plaza 
La  Jolla,  California  92037 

2  Systems,  Science  and  Software, 

Inc. 

ATTN:  Dr.  R.  Englemore 
Mr.  R.  Sklarew 
P.  0.  Box  1620 
La  Jolla,  California  92037 

1  TRW  Systems  Group 
ATTN:  Tech  Lib 
One  Space  Park 
Redondo  Beach,  California 
90278 

1  United  Aircraft  Corporation 
Research  Laboratory 
ATTN:  Mr.  R.  Bull is 
400  Main  Street 
East  Hartford,  Connecticut 
06108 

1  The  American  University 
Department  of  Physics 
ATTN:  Dr.  R.  Kay 
Massachusetts  and  Nebraska 
Avenues,  N.W. 

Washington,  D.C.  20016 

1  Cornell  Aeronautical  Lab.,  Inc. 
ATTN:  Mr.  R.  Fluegge 
P,  0,  Box  235 
Buffalo,  New  York  14221 


1  Georgia  Institute  of  Technology 
Chemistry  Department 
ATTN:  Dr.  T.  Moran 
Atlanta,  Georgia  30332 

1  Georgia  Institute  of  Technology 
Phys i cs  Department 

ATTN:  Dr.  E.  McDaniel 

Atlanta,  Georgia  30332 

2  The  Johns  Hopkins  University 
Chemistry  Department 

ATTN:  Dr.  W.  Koski 

Dr.  J.  Kaufman 
34th  and  Charles  Streets 
Baltimore,  Maryland  21218 

1  Kansas  State  University 
Physics  Department 
ATTN:  Prof.  L.  Leachmann 
Manhattan,  Kansas  66502 

1  Lincoln  Laboratory  (MIT) 

ATTN:  Mr.  J.  Panne 11 
244  Wood  Street 

Lexington,  Massachusetts  02173 

1  Lycoming  College 
Department  of  Physics 
ATTN:  Mr.  M.  Fineman 
Williamsport ,  Pennsylvania 
17701 

2  Midwest  Research  Institute 
ATTN:  Dr.  T.  Milne 

Dr,  F,  Greene 
425  Volker  Boulevard 
Kansas  City,  Missouri  64110 

1  State  University  of  New  York 
Department  of  Atmospheric 
Sciences 

ATTN:  Dr.  V.  Mohnen 
Albany,  New  York  12203 


66 


DISTRIBUTION  LIST 


No.  of 
Copies 


Organization 


No.  of 

Copies  Organization 


1  State  University  of  New  York 
at  Buffalo 

Department  of  Physics 
ATTN :  Mr .  G .  Brink 
Buffalo,  New  York  14214 

3  New  York  University 
Department  of  Physics 
ATTN:  Dr.  B.  Bederson 

Dr.  B.  Lippmann 
Dr.  T.  Miller 
University  Heights 
Bronx,  New  York  10453 

4  Pennsylvania  State  University 
Ionosphere  Research  Laboratory 
ATTN:  Prof.  A.  Waynick 

Dr.  J.  Gibbons 
Dr.  L.  Hale 
Dr.  M.  Zabielski 
University  Park,  Pennsylvania 
16802 

1  Princeton  University 
Department  of  Chemistry 
ATTN:  Prof.  M.  Kaufman 
Princeton,  New  Jersey  08540 

1  San  Die^o  State  College 
Department  of  Physics 
ATTN:  Mr.  F.  Wolf 
San  Diego,  California  92115 

1  Stanford  Research  Institute 
ATTN:  Mr.  J.  Norman 
4810  Bradford  Boulevard,  N.W. 
Huntsville,  Alabama  35805 

4  Stanford  Research  Institute 
ATTN:  Dr.  A.  Peterson 
Dr.  G.  Cook 
Dr.  D.  Lorents 
Dr.  W.  Aberth 
333  Revenswood  Avenue 
Menlo  Park,  California  94025 


4  Stanford  Research  Institute 


ATTN:  Dr. 

J. 

Peterson 

Mr. 

J. 

Lomax 

Dr. 

J. 

Moseley 

Dr. 

F. 

Smith 

333  Ravenswood  Avenue 
Menlo  Park,  California  94025 

1  University  of  Alaska 
Geophysical  Institute 
ATTN:  Tech  Lib 
College,  Alaska  99735 

1  University  of  California 
Department  of  Physics 
ATTN:  Prof.  L.  Loeb 
Berkeley,  California  94720 

1  University  of  California 
Department  of  Applied 
Electrophysics 
ATTN:  Mr.  P.  Banks 
P.  0.  Box  109 

La  Jolla,  California  92037 

6  University  of  Colorado 

Joint  Institute  for  Laboratory 
Astrophysics 
ATTN:  Dr.  Gordon  Dunn 
Dr.  S.  J.  Smith 
Dr.  A.  Phelps 
Dr.  G.  Sinnott 
Dr.  W.  Lineberger 
Dr.  L.  Kieffer 
Boulder,  Colorado  80304 

1  University  of  Colorado 
Laboratory  for  Atmospheric 
and  Space  Physics 
ATTN:  Mr.  J.  Pearce 
Boulder,  Colorado  80304 


67 


DISTRIBUTION  LIST 


No.  of 

Copies  Organization 


No.  of 

Copies  Organization 


2  University  of  Illinois 
Electrical  Engineering 

Department  Aeronomy  Laboratory 
ATTN:  Prof.  S.  Bowhill 
Prof.  C.  Sechrist 
Urbana,  Illinois  61803 

2  University  of  Maryland 

Institute  of  Molecular  Physics 
ATTN:  Prof.  J.  Vanderslice 
Prof.  T.  Wilkerson 
College  Park,  Maryland  20742 

1  University  of  Maryland 

Institute  of  Fluid  Dynamics 
and  Applied  Mathematics 
ATTN:  Dr.  T.  Rosenberg 
College  Park,  Maryland  20742 

6  University  of  Pittsburgh 
Cathedral  of  Learning 
ATTN :  Dr.  M.  Biondi 
Dr .  T .  Donahue 
Prof.  W.  Fite 
Prof.  E.  Gerjuoy 
Dr.  F.  Kaufman 
Dr.  E.  Zipf 
400  Bellefield  Avenue 
Pittsburgh,  Pennsylvania 
15213 


1  University  of  Texas 
Physics  Department 
ATTN:  Dr.  W.  Robertson 
Austin,  Texas  78712 

3  University  of  Texas 

Southwest  Center  for  Advanced 
Studies 

ATTN:  Dr.  F.  Johnson 

Dr.  C.  Collins 
Dr.  W.  Hurt 
P.  0.  Box  30365 
Dallas,  Texas  75230 

2  Wayne  State  University 
ATTN:  Dr.  A.  Hockstim 

Dr.  R.  H.  Kummler 
Department  of  Engineering 
Detroit,  Michigan  48202 

1  Yale  University 

ATTN:  Dr.  G.  J.  Schulz 

400  Temple  Street 

New  Haven,  Connecticut  06520 

Aberdeen  Proving  Ground 
Ch,  Tech  Lib 
Marine  Corps  Ln  Ofc 
CDC  Ln  Ofc 

CDCMA ,  ATTN:  Mr.  G.  Nelson 


68 


UNCLASSIFIED 

Security  Classification 


DOCUMENT  CONTROL  DATA  -  R  &  D 


>2«.  REPORT  SECURITY  CLASSIFICATION 


_ ( Security  ciaasi  iication  oi  title,  body  oi  abstract  and  indexing  annotation  must  be  entered  when  the  overall  report  ia  classified) 


1.  ORIGINATING  ACTIVITY  (Corporate  author) 

US  Army  Ballistic  Research  Laboratories 
Aberdeen  Proving  Ground,  Maryland  21005 


Unclassified 


NUMERICAL  MODEL  FOR  DRIFT  TUBE  ANALYSIS 


A.  DESCRIPTIVE  NOTES  (Type  of  report  and  inclueive  datee) 


G.  E.  Keller,  M.  R.  Sullivan,  L.  M.  Co lonna- Romano ,  and  M.  D.  Kregel 


e  REPORT  DATE 


•.  CONTRACT  OR  GRANT  NO. 


7 a.  TOTAL  NO.  OF  PAGES  |7b.  NO.  OF  REFS 


b  PROJEC 


tnoRDT5E  1T061102B53A/15 


REPORT  NO.  1617 


Ob.  OTHER  REPORT  NO(S)  (Any  other  numbers  that  may  be  aaalgyted 
this  report) 


10.  DISTRIBUTION  STATEMENT 


Distribution  limited  to  US  Government  agencies  only.  Other  requests  for  this 
document  must  be  referred  to  Director,  USA  Ballistic  Research  Laboratories,  ATTN: 

AMY  RR  _  Y^F  P  vr\  t  -  i  n  rr  (?Y»Atmrl  Uov.tr!  11  fine 

*  - - -  4  ~  a  V-*  w  *  *  aw  v  uiuuuu,  1'iai)  ianu  ^.xuuu. 


13.  ABSTRACT 


A  FORTRAN  IV  program  has  been  written  which  uses  an  iterative  technique  to  model  ion 
transDort  and  chemistrv  in  a  drift  tube.  The  urogram  miiitthatf  ..  =  »=  finite 

|  difference  calculus  to  generate  time  arrival  histories  at  several  points  in  the  drift 
|  space  for  the  case  of  three  ion  species  undergoing  populating  and  depopulating 

reactions  and  at  the  same  time  drifting  with  the  electric  field  and  diffusing  in  the 
radial  and  axial  directions.  To  facilitate  comparison  with  drift  tube  data,  the 


r nmm ifpr  nntnnt  inplnrlac  r.  1  nf c 

W*  vvr  A  VV*  X  IIV  XL4V*V^J  1^10  tJ 


1^  X  W11.L  iwLA  IICIIID  VCXZJUD  exilic. 


roan  i  a  -tr* 
UU  laviiR/J 


REPLACES  DD  FORM  1479.  I  JAN  E4.  WHICH  IE 
OBSOLETE  FOR  ARMY  USE. 


UNCLASSIFIED 

ecurity  Classification 


UNCLASSIFIED 

Security  Classification 


Ionic  Mobilities 
Diffusion  Coefficients 
Ion-Molecule  Reactions 
Drift  Tubes 
Drift  Tube  Models 
Drift  Tube  Analysis 


UINLLASbihihD 


Security  Classification 


